TECm'"-^ 

RE, 


library 

COPY 


5  0644  00004900  1 


m&2LL SjL 


TECHNIC 


AL  REPORT  NO.  11808 


VALUATION  OF  ACCURACY  OF  MEDIAN  RANKS  AND  MEAN 
RANKS  PLOTTING  FOR  RELIABILITY  ESTIMATION 
USING  THE  WEIBULL  DISTRIBUTION 


/ 


■  ^ 

i  Atoozotoioi'j 

-  by.  S.  B.  Catalano 

Materials  Branch 

TACOM 

VEHICULAR  COMPONENTS  &  MATERIALS  LABORATORY 


U.S.  ARMY  TANK  AUTOMOTIVE  COMMAND  Warren,  Michigan 


Distribution  of  this  document 
is  unlimited . 


V 

X 

CV) 

X? 


L 


The  citation  of  commercial  products 
in  this  report  does  not  constitute 
an  official  indorsement  or  approval 
of  such  products. 


The  findings  In  this  report  are  not  to  be 
construed  as  an  official  Department  of  the 
Army  position,  unless  so  designated  by  other 
authorized  documents. 


TECHNICAL  REPORT  NO.  11808 


EVALUATION  OF  ACCURACY  OF  MEDIAN  RANKS  AND 
MEAN  RANKS  PLOTTING  FOR  RELIABILITY  ESTIMATION 
USING  THE  WEIBULL  DISTRIBUTION 

By 

S.  B.  Catalano 

June  1973 


AMCMS  CODE:  611101.11.8UU00 


MATERIALS  BRANCH 


Reproduced  From 
Best  Available  Copy 


ABSTRACT 


Estimates  of  the  Weibull  distribution  parameters  were  made  employing  the 
mean  ranks  estimator;  the  estimates  were  repeated  using  the  median  ranks 
estimator.  These  estimates  were  compared  to  known  values  of  the  Weibull 
distribution  parameters.  This  made  it  possible  to  compare  the  results 
obtained  using  either  estimator  (mean  ranks  or  median  ranks)  and  to 
determine  the  relative  merits  of  using  either  estimator.  This  study  made 
use  of  a  digital  computer  and  employed  Monte-Carlo  techniques  to  simulate 
Weibull  distributed  failure  times.  These  failure  times  may  represent 
tank-automotive  component  failures. 
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INTRODUCTION 


One  of  the  earliest  applications  of  the  Weibull  distribution  in  this 
country  was  in  1951  in  a  paper  presented  by  Professor  Weibull  (l). 

Its  use  since  then  has  been  predominantly  in  the  analysis  of  life  test 
data  in  which  the  variable  of  interest  is  lifetime,  t.  The  Weibull 
distribution  is  of  interest  to  TACOM  from  the  standpoint  of  analysis  of 
life  test  data  on  tank- automotive  components.  The  Weibull  density 
function  for  the  random  variable,  t,  is: 

f (t)  =  1  exp 

f(t)  =  0  for 

where  ft  is  the  shape  parameter,  or  slope  parameter  (usually  the  value  of 
ft  is  near  3.5  for  tank-automotive  components;  this  means  that  the  failure 
distribution  curve  has  the  familiar  bell  shape). 

■&  is  the  scale  parameter,  or  characteristic  life  parameter  (the 
units  of  ■&  for  tank-automotive  components  are  usually  measured  in  miles 
or  cycles  till  failure). 

cK  is  the  location  parameter,  or  minimum  life  parameter. 

In  general  usage,  cK  =  0,  in  which  case: 

f(t)  =  /L  1  exp 

f(t)  =  0  for  t<0. 


In  this  study  we  will  be  concerned  with  the  cumulative  Weibull  distribution 
which  is  expressed  mathematically  as: 


F(t) 


or,  if  cK-  0: 


F(t)  =  1  -  exp 
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In  the  application  of  the  Weibull  distribution  in  the  analysis  of  life 
test  data,  use  is  made  of  estimators  for  the  value  of  F(t).  Estimates 
are  necessary  since  life  testing  data  yields  only  values  of  failure  time 
t;  the  values  of  f!>  and  -&-in  equation  (l)  remain  unknown.  Two  different 
estimators  used  for  this  purpose  are  called  "median  ranks"  and  "mean 
ranks"  estimators. 

Background  information  is  presented  in  Appendix  I.  The  computer  program 
flow  chart  is  shown  in  Figure  1.  Figures  2  through  11  in  Appendix  I 
present  various  aspects  of  the  Weibull  distribution. 


OBJECT 


The  objective  of  this  study  is  to  compare  the  accuracy  of  the  two 
commonly  used  estimators  called  "median  ranks"  and  "mean  ranks"  when 
employed  in  Weibull  distribution  failure  analysis  for  estimating  the 
Weibull  oarameters  and  -fr  . 


SUMMARY 


Weibull  distributed  failure  times  were  simulated  on  a  computer  via 
Monte-Carlo  techniques.  Using  these  simulated  values  of  failure  times, 
estimates  of  the  Weibull  parameters^  and  &■  were  computer  calculated 
using  the  median  ranks  estimator;  the  calculations  were  repeated  using 
the  mean^anks  estimator.  The  resulting  estimates  of  $  and  -e-vere 
labeled^1  and  #-and  were  compared  to  the  known  values’ of  £  and  £-that 
were  used  to  simulate  the  Weibull  distributed  failures.  In  this  manner, 
it  was  possible  to  compare  the  accuracy  of  median  ranks  and  mean  ranks 
estimators  in  calculating  estimated  values  of  Weibull  parameters  £  and  &  . 
The  study  was  carried  out  for  various  sample  sizes  and  for  various  values 
of  the  parameters  f5  and  & ;  various  degrees  of  suspended  data  were 
employed. 


CONCLUSIONS 


It  is  concluded  from  this  study  that ^se  of  the  median  ranks  approximator 
provides  a  better  approximation  for  @  than  mean  ranks  over  the  range  of 
parameters  employed  in^his  study.  This  also  tends  to  be  the  case  for 
estimating  the  slope,  {}  ,  for  cases  where  the  number  of  test  samples  is 
high;  this  tendency  is  enhanced  as  the  degree  of  suspended  data  is 
decreased. 
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The  mean  ranks  estimator  provides  a  better  estimate  of  p  than  the  median 
ranks  estimator  only  in  the  region  of  small  sample  size  and  high  degree 
of  suspended  data.  But  this  region  is  a  region  wh^ch  is  inherently 
plagued  with  a  high  degree  of  error  in  estimating  regardless  of  which 
estimator  is  employed.  The  additional  error  that  is  encountered  by  using 
median  ranks  rather  than  mean  ranks  in  this  region  is  small  in  comparison 
with  the  error  inherently  encountered  in  this  region. 


RECOMMENDATIONS 


It  is  deemed  advantageous  to  routinely  employ  the  median  approximator 
when  estimating  the  Weibull  parameters  @  and  ■&-  .  The  slight  loss  of 
advantage  in  the  region  of  small  sample  size  and  large  degree  of  suspended 
data  in  the  case  of  estimating  ft  can  either  be  corrected  for  by  reference 
to  Figures  llj  through  19  of  Appendix  IV  or  dismissed  as  being  negligible. 


TEST  PROCEDURE 


Theintent  of  this  study  was  to  objectively  compare  the  resulting  values 
of  ft  and  £  when  using  mean  ranks  with  those  values  of  ft  and  ^  obtained 
when  using  median  ranks,  for  both  suspended  data  tests  and  non-suspended 
data  tests.  The  comparison  was  made  by  measuring  ft  and  ^  in  each  of 
twenty  separate  simulated  failure  tests  using  first  the  mean  ranks 
estimator,  then  using  the  median  ranks  estimator.  The  standard  error 
statistic  was  employed  as  the  measure  of  accuracy  for  comparison  purposes. 
In  this  case,  the  mathematical  expressions  for  standard  error  in  the 
measurement  of  R  and  -0-  are  respectively 
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This  study  was  carried  out  for  a  wide  range  of  values  of  the  parameters 
@ ,  -£-and  sample  size,  n,  under  suspended  and  non-suspended  data  testing. 
The  values  of ft  used  were  0.9,  1,  2,  U  and  8.  The  values  of  &- used  were 
9,000,  10,000,  20,000,  U0,000,  80,000  and  160,000  miles.  Test  sample 
sizes  used  were  n  =  9j  10,  20,  I4O,  80  and  160  samples.  Tests  were  run 
for  60$,  h0%  and  0 %  (i.e.  non-suspended)  suspended  data  runs. 


The  entire  study  was  a  computer  study.  Failures  were  simulated  on  a 
digital  computer  using  Monte-Carlo  techniques.  Rather  than  plotting  the 
resulting  points  and  drawing  the  best  fit  straight  line  through  these 
data  points  by  sight,  a  computer  subroutine  was  used  to  compute  the  best 
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fit  using  least  squares  fitting.  In  this  manner  it  is  expected  that  the 
results  obtained  are  objective  and  free  of  differences  due  to  personal 
traits  in  line  plotting  by  eye.  The  computer  was  then  asked  to.cornpute 
0  and  ^and  then  to  compute  and  print  the  standard  error  [_2(  P  ^ 

and  £-)2/20 7 1/2  for  each  of  the  three  degrees  of  suspended  data 

tests  employed  ancl  for  each  combination  of  parameters  0  ,  and  n  employed 
in  this  study. 

The  flow  diagram  of  the  computer  program  used  in  this  study  is  shown  in 
Figure  1.  Each  step  in  the  flow  diagram  will  be  explained  in  detail. 

The  computer  program  as  listed  just  prior  to  computer  execution  is  shown 
in  Appendix  II.  The  resulting  printout  is  shown  in  Appendix  III. 

Weibull  distributed  failure  times  are  simulated  (via  Monte-Carlo  method) 
by  evaluating  equation  (l)  for  t  using  random  numbers  from  0.0  to  1.0 
for  F(t).  In  real  life,  the  selection  of  random  numbers  from  0.0  to  1.0 
represents  the  random  selection  of  test  samples  from  the  entire  population 
of  samples.  The  significance  of  the  random  numbers  ranging  from  0.0  to 
1.0  is  that  the  test  samples  selected  may  equally  well  be  the  first  to 
fail,  the  last  to  fail  or  to  fail  at  any  time  between  the  first  and  last 
failures.  In  other  words,  the  values  from  0.0  to  1.0  represent  a  ranking 
of  the  failures  onto  a  percentage  scale.  In  equation  (l),  F(t)  represents 
this  ranking.  Evaluating  equation  (l)  for  t  using  the  random  numbers  for 
values  of  F(t)  has  the  effect  of  grouping  or  modulating  the  failure  times 
such  that  they  will  be  Weibull  distributed  random  failure  times.  In 
other  words,  the  failure  times  so  generated  simulate  failure  times  that 
would  occur  for  samples  that  fail  according  to  the  Weibull  probability 
distribution. 

The  need  for  random  numbers  in  this  study  is  apparent.  Random  number 
generators  are  not  available  in  all  computers.  The  computer  used  in  this 
study  was  a  Control  Data  Corporation  computer  CDC  Model  6600.  It 
possesses  a  system  function  RAKF(O)  which  provides  a  random  number  from 
0.0  to  1.0  each  time  RANF(O)  is  requested  by  the  program.  The  number  of 
decimal  places  utilized  in  this  study  for  RANF(O)  was  eight  places. 

In  this  study  sets  of  random  numbers  are  used.  The  CDC  system  function 
RANF(O)  furnishes  a  new  set  of  random  numbers  each  time  a  new  set  is 
called  for  in  a  program.  However,  each  time  the  program  is  resubmitted 
to  the  computer,  it  gives  the  exact  same  sets  of  numbers  in  exactly  the 
same  sequence  from  start  to  finish. 

The  first  step  in  the  program  flow  diagram  indicates  the  calling  for  a 
set  of  n  random  numbers. 
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STEP 

NUMBER 

1. 

2. 

3. 

U. 

6. 

7. 

8. 

9. 

10. 

11. 

12. 

13. 

1U. 

15. 

16. 

17. 

18. 


OPERATION 

Gall  for  set  of  n  random  numbers 
Rank  order  the  set  of  n  random  numbers 
Generate  n  Monte-Carlo  Simulated  Failure  Times 
Generate  n  Median  ranks  (rank  ordered) 

Generate  n  Mean  ranks  (rank  ordered) 

Perform  transformation  of  failure  times  to  Weibull  Probability  Axes 

Perform  transformation  of  mean  ranks  to  Weibull  Probability  Axes 

Perform  transformation  of  median  ranks  to  Weibull  Probability  Axes 

Pair  the  n  transformed  failure  times  with  the  n  transformed  mean 
ranks 

Pair  the  n  transformed  failure  times  with  the  n  transformed 
median  ranks 

^  y\ 

Determine  p  and  &  for  mean  ranks  using  least  squares  fitting 
subroutine 

^  A 

Determine  p  and  for  median  ranks  using  least  squares  fitting 
subroutine 

Repeat  steps  1  thru  12  M  times  (where  M  =  20) 

Determine  standard  errors  jjE(/^  &  jjl(£-  0-)^/mP^ 

for  mean  ranks  L“  L.  -J 

Determine  standard  errors 

for  median  ranks  u  '  *“ 

Change  to  next  value  of  n  and  repeat  steps  1-15 

Change  to  next  value  of  0  and  repeat  steps  1-16 

Change  to  next  value  of  -0-  and  repeat  steps  1-17 


FIGURE  1.  Computer  Program  Flow  Chart 
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The  second  step  in  the  program  flow  diagram  shows  that  the  set  of  n  random 
numbers  generated  in  the  previous  step  will  be  rank  ordered.  The  ordering 
is  from  lowest  number  to  highest  number.  The  random  numbers  are  then  used 
in  this  order  to  generate  the  simulated  failure  times.  The  resulting 
failure  times  are  thus  generated  in  rank  (numerical  or  chronological) 
order.  Aside  from  it  being  psychologically  satisfying  to  have  the  failure 
times  occurring  in  chronological  order,  thereby  giving  the  simulation  a 
real-life  flavor,  it  facilitates  pairing  of  failure  times  with  the  proper 
median/mean  ranks:  the  first  failure  time  with  the  first  value  of 
median/mean  ranks;  the  second  failure  time  with  the  second  value  of 
median/mean  ranks,  etc.,  till  the  last  failure  time  is  paired  with  the 
last  value  of  median/mean  ranks.  This  pairing  is  for  purposes  of  forming 
coordinate  points  to  be  plotted/fitted  onto  Weibull  probability 
coordinate  axes.  Note  that  it  is  immaterial  whether  the  rank  ordering 
step  is  performed  before  or  after  the  generation  of  failure  times  or  for 
that  matter  whether  it's  done  at  all.  What  is  important  is  that  somehow 
the  jth  failure  time  is  paired  with  the  median/mean  ranks.  Rank 
ordering  facilitates  this  pairing. 

The  third,  fourth  and  fifth  steps  in  the  program  flow  diagram  indicate 
the  generation  of  a  set  of  n  rank  ordered  failure  times,  median  ranks  and 
mean  ranks  respectively. 

Up  to  this  point,  the  computer  has  generated  a  set  of  n  rank  ordered 
random  numbers,  a  set  of  n  rank  ordered  failure  times,  a  set  of  n  rank 
ordered  median  ranks  and  a  set  of  n  rank  ordered  mean  ranks.  The  next 
step  is  to  perform  a  transformation  of  axes  on  the  above  sets  of  numbers 
to  Weibull  probability  axes,  (X,  Y).  The  failure  times,  t,  are  transformed 
according  to  equation  (6)  of  Appendix  I.  Since  the  random  numbers 
represent  various  values  of  F(t)  and  the  median  ranks  and  the  mean  ranks 
are  estimates  of  F(t),  these  three  sets  are  transformed  according  to 
equation  (5>)  of  Appendix  I.  Steps  6,  7  and  8  of  the  program  flow  diagram 
indicate  these  transformations  of  axes. 

The  next  block  of  steps  in  the  program  flow  diagram  is  for  performing  the 
least  squares  fitting  of  a  straight  line  to  the  set  of  coordinate  points 
on  Weibull  probability  coordinate  axes.  The  first  step  in  this  process 
is  to  pair  the  jth  transformed  failure  times  with:  l)  the  median 
ranks,  and  2)  with  the  jth  mean  ranks;  this  is  shown  in  steps  9  and  10. 

The  next  step  is  to  perform  the  actual  least  squares  fitting.  This  is 
done  with  the  subroutine  LSFIT.  It  makes  use  of  the  least  squares 
fitting  procedure  as  outlined  by  Cullity  (2).  Accordingly,  subroutine 
LSFIT  computes  the  first  and  second  normal  equations  for  the  simulated 
points  on  the  Weibull  probability  coordinate  axes.  These  two  equations 
are  linear  equations;  simultaneous  solution  yields  the  slope  of  the 
straight  line  and  the  Y-intercept.  In  this  study,  the  slope  is^  and  the 
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Y-intercept  is  -p\n-&  from  which  the  value  of  -&  is  readily  obtained. 
Solution  of  the  simultaneous  equations  is  done  in  the  main  program.  ^ 
Steps  11  and  1?  of  the  program  flow  diagram  indicate  computation  of  p  and 
for  the  median  and  mean  ranks. 


The  next  step  in  the  program  flow  diagram  indicates  that  all  of  the 
previous  steps  are  to  be  repeated  20  times.  In  real  life,  this  would 
simulate  repeating  20  times  the  entire  process  of  randomly  selecting  n 
test  samples  to  be  life  tested,  plotting  the  failure  times  (on  Weibull 
probability  paper)  vs.:  l)  mean  ranks,  and  2)  median  ranks,  and  then 
drawing  the  best  fit  straight  line  for  both  cases  and  determining  p  and  ^ 
from  the  slopes  and  intercepts  of  these  lines.  This  results  in  20  values 
of  and  ^for  mean  ranks  case  and  20  for  the  median  ranks  case. 


In  the  next  step  shown  in  the  flow  diagram,  the  program  computes  the 
standard  errors  )^/2CQ  and  -  &')w20j  V2  from  the  20 

values  of  fi  and  ^arrived  at  in  the  previous  step.  This  is  done  for  both 
the  median  ranks  and  mean  ranks  values  of  ft  and  & ,  and  for  each  of  the 
three  degrees  of  suspended  data  (60$,  k0%  and  0%).  This  calculation  is 
performed  by  the  standard  error  subroutine  MDEV  listed  at  the  end  of  the 
program.  / 

As  a  reminder,  it  should  be  noted  here  that  the  values  of  and  #-are 
known.  They  were  used  early  in  the  program  in  equation  (l)  to  generate 
the  simulated  failure  times. 


The  next  three  steps  in  the  program  direct  the  computer  to  repeat  the 
computation  for  new  parameter  values,  changing  the  value  of  one  parameter 
at  a  time  until  all  possible  combinations  of  parameter  values  ofy#  ,  and 
n  used  in  this  study  have  been  used. 


RESULTS  AND  DISCUSSION 

The  numerical  values  of  standard  errors  and 

DECS--e-)2/2qjV  ^  as  calculated  and  printed  by  the  computer  are  shown 
in  Appendix  III.  These  values  of  standard  error  are  tabulated  according 
to  degree  of  suspended  data  (i.e.  60%,  U0%  or  0%)  and  are  grouped  according 
to  parameter  values  of and-#-  ,  each  group  making  use  of  each  value  of 
n  (sample  size)  employed  in  this  study.  There  are  2,160  entries  of  standard 
error  charted,  thus  making  it  overwhelmingly  difficult  to  interpret 
simply  by  comparing  numerical  values.  For  this  reason,  these  values  of 
standard  error  were  plotted;  they  were  plotted  vs.  sample  size,  n.  These 
plots  are  shown  in  Appendix  IV  (Figures  lit  thru  25>) .  Each  figure  contains 
plots  for  60%,  U0%  and  0%  suspended  data  thus  facilitating  interpretation 
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of  results  as  a  function  of  degree  of  suspended  data.  These  plots  are 
drawn  for  each  of  the  five  values  of  $  used  in  this  study  and  for  both 
estimators  (median  ranks  and  mean  ranks).  In  this  manner,  the  results 
can  be  more  readily  interpreted  as  a  function  of  as  well  as  a  function 
of  estimator  employed.  Only  one  value  of  -&■  is  used  per  figure.  Inter¬ 
pretations  as  a  function  of  -Q-  must  therefore  be  made  by  making  comparisons 
between  the  various  figures. 

In  the  cases  of  60$  and  1*0$  suspended  data,  the  plots  do  not  include  the 
points  for  n  =  5  since  in  these  cases,  it  would  mean  that  the  tests  were 
terminated  after  only  two  and  three  samples,  respectively,  had  failed, 
which  is  much  too  small  a  number  of  failures  from  which  to  obtain  meaningful 
results.  Nevertheless,  the  computer  was  asked  to  perform  these  calculations. 
In  most  cases,  the  results  so  obtained  were  either  off  the  scale  of  the 
graphs  shown  or  larger  than  the  largest  number  allowed  to  be  printed  by 
the  computer  program  format  (in  this  case,  an  asterisk  is  shown  in  the 
computer  printout),  thereby  preventing  plotting  of  such  points. 

Figures  lli  through  19  of  Appendix  IV  are  plots  of  standard  error  encountered 
in  calculating^  ;  figures  20  through  25  of  Appendix  IV  are  for  standard 
error  encountered  in  calculating  . 

In  general,  the  two  effects  that  these  graphs  make  most  apparent  are  effects 
that  are  intuitively  expected;  these  are  that: 

1.  Standard  error  decreases  with  increasing  sample  size.  / 

2.  Standard  error  increases  with  increasing  degree  (or  percentage) 
of  suspended  data. 

Another  effect  that  is  immediately  apparent  frcrni  these  graphs  in  the  case 
of  standard  errors  encountered  in  calculating  (see  Figures  lli  thru  19) 
is  that  standard  error  increases  as  increases.  It  is  also  seen  in  this 
case  that  has  little  or  no  effect  on  standard  error.  However,  in  the 
case  of  standard  errors  encountered  in  calculating (see  Figures  20  thru 
25),  standard  error  decreases  as  ft  increases,  and  is  significantly 
increased  as-O-is  increased. 

The  zig-zagged  appearance  of  the  plots  is  due  to  connecting  the  plotted 
points  with  straight  lines.  Generally,  the  plots  exhibit  decreasing 
standard  error  with  increasing  sample  size  n,  however,  occasionally  a 
point  will  deviate  from  this  trend  and  will  be  higher  than  the  previous 
point,  when  it  was  expected  to  be  lower  than  the  previous  point.  This 
possibly  is  caused  by  random  sampling;  if  so,  it  is  a  reflection  of  the 
fluctuations  experienced  in  real  life  failure  testing  when  test  samples 
are  randomly  chosen  from  the  entire  nonulation  of  samples.  This  real  life 
"flavor"  is  reflected  into  the  results  of  this  simulation  study  since  a 
random  number  generator  was  used  in  generating  Weibull  distributed  failure 
times. 
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Solid  lines  in  the  plots  are  used  to  illustrate  the  results  obtained  when 
using  the  median  ranks  estimator;  broken  lines  are  used  for  the  results 
obtained  when  using  the  mean  ranks  estimator. 

The  following  observations  are  made  (relative  to  the  solid  vs.  broken 
lines)  in  the  case  of  standard  errors  encountered  in  calculating  ft  (refer 
to  Figures  lh  thru  19): 

1.  Sometimes  the  solid  lines  are  above  the  broken  lines;  sometimes 
they  are  below.  Sometimes  the  two  are  overlapping. 

2.  The  broken  lines  tend  to  be  lower  than  the  solid  lines  for  lower 
values  of  n,  and  tend  to  be  higher  for  higher  values  of  n. 

3.  The  solid  lines  tend  to  be  lower  than  the  broken  lines  as 
increases  and  as  the  degree  (or  percentage)  of  suspended  data  is 
decreased. 

h.  The  solid  lines  cross  or  touch  the  broken  lines  at  least  once  in 
the  range  of  values  of  n  used  in  this  study. 

The  following  observations  are  made  (relative  to  the  solid  vs.  broken 
lines)  in  the  case  of  standard  errors  encountered  in  calculating  4^( refer 
to  Figures  20  thru  2$) : 

1.  The  solid  lines  are  either  below  the  broken  lines  or  the  two  are 
overlapping.  In  most  cases,  even  when  the  two  lines  are  shown 
overlapping,  reference  to  the  numerical  values  of  the  plotted 
points  reveals  that  the  solid  line  has  the  lower  value  at  that 
point.  There  were  several  exceptions  to  this.  However,  in 
these  cases  the  two  values  were  so  close  numerically  that  it  was 
not  possible  to  depict  the  difference  when  plotting  these  values. 

2.  The  solid  and  broken  lines  tended  to  be  closer  together  as 

l)  increased,  2)  the  degree  of  suspended  data  decreased,  and 
3)  n  increased. 

As  pointed  out  in  Appendix  I  of  this  report,  L.  G.  Johnson  prefers  using 
median  ranks  rather  than  using  mean  ranks  when  drawing  the  best  fit  line 
through  the  data  points  on  Weibull  probability  paper  by  eye.  He  points 
out  that  by  using  median  ranks  one  avoids  underestimating  the  slope,  ft  , 
of  the  line  in  the  region  of  the  lower  extreme  of  the  graph  (i.e.  in  the 
region  where  only  relatively  few  of  the  test  samples  have  failed).  In 
terms  of  degree  of  suspended  data,  this  corresponds  to  the  region  of 
higher  degree  (or  percentage)  of  suspended  data. 


The  results  of  this ^tudy  indicate  that  the  tendency  to  err  in  the 
estimate  of  slope,  ,  at  the  higher  degrees  of  suspended  data  tends  to 
vary  with  n,  the  test  sample  size.  The  tendency  is  that  for  larger 
'sample  sizes,  median  ranks  are  more  accurate  than  mean  ranks,  whereas 
for  smaller  sample  sizes,  the  reverse  is  true. 

Although  Mr.  Johnson  did  not  expressly  discuss  the  relative  merits  of 
using  median  ranks  vs.  mean  ranks  when  estimating^,  one  must  assume 
that  he  prefers  median  ranks  in  this  case  also,  since  an  error  in  slope, 
p  ,  would  mathematically  reflect  an  error  in  the  estimate  of  The 

results  in  this  respect  in  the  present  study  indicate  that  use  of  median 
ranks  is  as  accurate  as,  or  more  accurate,  than  use  of  mean  ranks. 


\ 
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A  short  discussion  of  the  general  properties  of  the  Weibull  distribution 
follows : 


W 


The  expected  value  of  the  Weibull  distribution  is 


E(t)  =-0-p 


The  variance  of  the  Weibull  distribution  is 


V(t) 


(2) 


(3) 


The  shape  of  the  Weibull  density  function 

"changes  from  a  highly  positively  skewed  distribution 
when  0  =  0.5,  to  a  simple  exponential  distribution 
when  0  =  1,  to  an  essentially  Gaussian  normal  distribution, 
when  0  -  3.5,  to  a  negatively  skewed  distribution  when 
0  =  6  or  more. "  (3) 

These  various  shapes  are  shown  in  Figure  2. 

When  P)  -  1, 


f(t)  =  i  exp 

-& 

which  verifies  that  when  ft  -  1,  the  Weibull  distribution  is  equivalent 
to  the  simple  exponential  distribution.  Figure  3  illustrates  that  when 
0  -  3.5,  the  Weibull  distribution  is  essentially  a  Gaussian  normal 
distribution.  This  illustration  was  made  by  comparing  a  plot  of  the 
standard  normal  curve  with  a  plot  of  the  Weibull  distribution  having 
fi  =  3*5  and  variance  equal  unity.  Note  that <X =  0  in  this  plot  of  the 
Weibull  distribution,  and  that  the  plot  would  "slide"  to  a  new  position 
along  the  horizontal  axis  as  various  non-zero  values  of<7<are  used. 
Figure  3  also  illustrates  a  plot  of  the  Weibull  distribution  for0-  3*5 
with  variance  equal  U  (i.e.  standard  deviation  =2). 


t 
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Figures  h  through  8  are  plots  of  the  Weibull  distribution  that  were  made 
using  values  of  equal  to  3«0,  3*3>  3*5,  3»7  and  U.O.  In  each  case  the 
variance  equals  unity.  For  comparison  purposes  a  standard  normal  curve 
is  also  plotted  on  each  figure.  This  series  of  plots  was  made  to  verify 
that  the  Weibull  distribution  having^  =3*5  conforms  closest  to  the 
standard  normal  curve.  It  is  seen  from  these  figures  that  of  those  values 
of  fe  used,  the  Weibull  distribution  having^  =  3*5  conforms  closest  to 
the  standard  normal  curve.  The  greatest  discrepancy  here  between  the  two 
distributions  is  observed  to  be  at  the  peaks  of  the  distributions. 
Nevertheless,  it  is  apparent  that  the  Weibull  distribution  having  ft  =  3*5 
and<5"  =  1  is  a  good  approximation  to  the  standard  normal  distribution. 

Figures  9,  10  and  11  illustrate  that  the  plot  of  the  Weibull  distribution 
becomes  taller,  slimmer  and  more  negatively  skewed  as  0  gets  larger.  In 
each  of  these  plots,  the  variance  equals  unity;  the  values  of  in  Figures 
9,  10  and  11  are  respectively  10,  I4O  and  100. 

From  equations  (?)  and  (3),  it  is  seen  that  as 0  approaches  infinity,  the 
expected  value  equals  fr  and  that  Q~  equals  infinity.  It  appears  that  the 
Weibull  distribution  is  of  limited  use  at  the  larger  values  of  (such 
as  ft  =  100) .  Its  usefulness  lies  primarily  in  the  region  of  the  lower 
values  of  ft  ,  especially  at  fi  =1,  where  the  Weibull  distribution  is 
equivalent  to  the  exponential  distribution,  and  at  =  3.5  where  it  is 
essentially  equivalent  to  the  Gaussian  normal  distribution.  Both  the 
exponential  distribution  and  the  Gaussian  distribution  are  useful  in 
failure  analysis.  Consequently  the  Weibull  distribution,  due  to  its 
versatility,  is  very  useful  in  failure  analysis.  In  fact,  its  usefulness 
transcends  that  of  either  the  exponential  or  Gaussian  distribution  since 

"the  exponential  distribution  is  applicable  as  a  model  for 
failure  times  only  if  the  failure  rate  is  constant  over  time.  In 
reality,  failure  rates  which  change  with  time  are  sometimes 
encountered.  The  normal  distribution  is  a  realistic  model  only 
if  an  increasing  failure  rate  is  encountered.  The  Weibull 
distribution  is  continuous  and  can  account  for  a  decreasing 
failure  rate."  Qj) 

To  further  illustrate  the  usefulness  of  the  Weibull  distribution,  reference 
is  made  to  Figure  12.  It  illustrates  the  typical  life  history  of  a 
population  of  units  of  a  complex  product.  The  initial  or  "de-bugging" 
phase  is  caused  by  the  short  life  of  marginal  units;  this  phase  is 
characterized  by  a  high  but  decreasing  failure  rate.  This  phase  of  the 
life  history  can  be  handled  with  the  Weibull  distribution  using ^ ■Cl, 

The  second  phase  is  characterized  by  a  low  and  relatively  constant  failure 
rate  which  lasts  till  the  units  begin  to  wear  out.  The  failure  rate 
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FIGURE  12.  Typical  life  history  of  a  population  of  units  of  a  complex 
product  (not  to  scale  regarding  phase  lengths). 
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during  this  period  is  due  to  chance  failures.  This  portion  of  the  life 
history  can  be  handled  with  the  Weibull  distribution  using  =1.  The 
last  phase  is  the  wear  out  phase  and  is  characterized  by  an  increasing 
failure  rate.  This  portion  of  the  life  history  can  be  handled  with  the 
Weibull  distribution  using^>l.  (5)  It  is  seen  then  that  the  Weibull 
distribution  can  be  employed  in  failure  analysis  over  the  entire  life 
history  of  the  product,  and  is  consequently  of  considerable  value  in 
failure  analysis. 

The  computer  was  employed  in  this  study  to  simulate: 

1.  Weibull  distributed  failures  via  Monte-Carlo  techniques,  and 

2.  Use  of  Weibull  probability  paper  for  estimating  the  values  of 
P  and 

A  description  of  the  Monte-Carlo  technique  was  given  earlier  in  this 
report  as  part  of  the  section  on  Test  Procedures.  The  following  is  a 
description  of  the  use  of  Weibull  probability  paper  for  estimating 
and  . 

This  method  for  estimating  0  and  -^"employs  graph  paper  called  Weibull 
probability  paper.  The  graph  paper  has  its  axes  so  graduated  that  when 
Weibull  distributed  failures  are  plotted,  a  straight  line  plot  will 
result.  An  example  of  Weibull  probability  paper  is  shown  in  Figure  13. 
Failure  times,  t,  are  plotted  along  the  horizontal  axis.  The  vertical 
axis  represents  the  cumulative  Weibull  probability,  F(t)  (see  equation 

(D). 

Actual  life  testing  consists  of  selecting  n  test  samples  from  a  large 
but  unspecified  sized  population,  and  running  them  to  failure.  The  only 
data  that  is  obtained  from  such  a  test  is  the  various  values  of  time  to 
failure,  tj,  where  j  *  1  for  the  first  failure,  j  =  2  for  the  second 
failure,  etc.,  up  to  j  =  n  for  the  last  failure.  These  values  of  failure 
times  are  plotted  along  the  horizontal  axis.  The  corresponding  values 
of  F(tj)  are  obtained  from  an  estimator  and  are  plotted  along  the 
vertical  axis.  The  most  commonly  used  estimator  for  this  purpose  is 
called  the  "mean  ranks".  Algebraically,  it  is  expressed  as  j/(n  +  l). 
Using  the  symbol  above  a  quantity  to  indicate  "estimated  value  of", 
we  have 


F(tj)  =  j/(n  +  1) 


for  the  mean  ranks  estimation.  The  physical  meaning  of  "mean  ranks"  is 
that  if  the  same  failure  analysis  test  were  conducted  many  times,  the 
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mean  value  for  a  particular  F(tj)  would  be  j/(n+l).  That  is: 

the  mean  value  for  F(t^)  would  be  l/(n  +  l), 
the  mean  value  for  F(tp)  would  be  2/(n  +1), 


the  mean  value  for  F(tn)  would  be  n/(n  +  l). 

The  derivation  of  the  expression  j/(n  +  l)  for  the  mean  ranks  is  shown 
by  L.  G.  Johnson  (6)  and  L.  R.  Lamberson  (?)•  Messrs  Kao  and  Goode  (8), 
(9),  (10)  have  done  much  work  with  the  Weibull  distribution;  much  of 
their  studies  make  use  of  mean  ranks.  The  United  States  Army  (U) 
prescribes  the  use  of  mean  ranks  in  connection  with  life  testing  using 
the  Weibull  distribution. 

Another  estimator  used  for  this  purpose  is  called  the  "median  ranks". 

The  exact  values  of  the  median  ranks  are  obtained  from  tables,  however, 
an  approximating  equation  for  the  median  ranks  is  given  by  E.  J.  Gumbell 
(ll)  as  (j  -  .3)/(n  +  .It).  Using  it  we  have 

F(tj)  =  (j  -  .3)/(n  +  .lj) 

for  the  median  ranks  estimation.  The  physical  meaning  of  "median  ranks" 
is  that  if  the  same  failure  analysis  experiment  were  conducted  many  times, 
the  median  value  for  a  particular  F(tj)  would  be  (j  -  .3)/(n  +  .U). 

That  is 

the  mean  value  for  F(t^)  would  be  (l  -  .3)/(n  +  .U), 
the  mean  value  for  F(t2)  would  be  (2  -  .3)/(n  +  .lj), 


the  mean  value  for  F(tn)  would  be  (n  -  .3)/(n  +  .[(). 

Mr.  L.  G.  Johnson  has  written  many  papers  (6),  (12)  and  books  (13),  (lh) 
dealing  with  analysis  using  the  Weibull  probability  distribution.  He 
favors  use  of  median  ranks  when  hand  drawing  the  straight  line  on 
Weibull  probability  paper.  He  writes  (l5): 
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"To  draw  the  line  very  near  the  lower  extreme  values  when  mean 
ranks  are  used  would  lead  to  a  slope  which  is  too  small, 
because  of  the  high  probability  of  lower  extreme  values 
falling  considerably  to  the  left  in  such  a  case.  On  the 
other  hand,  it  has  been  found  that  if  median  ranks  are  used 
in  plotting,  then  the  danger  of  under  estimating  the  slope 
is  eliminated  because,  in  this  case,  a  point  is  just  as 
liable  to  fall  to  the  right  as  to  the  left  of  the  maximum 
likelihood  line.  In  other  words,  we  can  very  quickly  arrive 
at  an  estimate  of  the  population  parameters  by  drawing  a 
line  by  sight  which  takes  the  general  direction  of  the  array 
of  points  and  which  splits  the  array  5>0-£0.  For  this  reason 
median  ranks  are  preferred  to  mean  ranks." 

It  should  be  noted  that  regardless  of  which  estimator  (mean  ranks,  or 
median  ranks)  is  used,  0  F(tj)  —  1.  In  terms  of  percentages, 

F(tj)  lies  between  0 %  and  10052.  In  the  life  test  experiment  described 
above,  F( t j )  represents  percentage-wise,  the  number  of  samples  out  of 
the  entire  population  of  samples  that  will  have  a  lifetime  less  than  or 
equal  to  the  lifetime  of  the  test  sample  to  fail. 

If  it  is  not  possible  to  draw  a  straight  line  through  the  points  plotted 
on  Weibull  probability  paper,  then  either  the  failures  are  not  Weibull 
distributed  or  the  minimum  life  parameter,  cK. ,  is  not  zero.  If  the 
failures  are  not  Weibull  distributed,  then  Weibull  probability  paper 
cannot  be  used.  If  o(  is  not  zero,  the  data  can  be  adjusted  via  a 
transformation  of  the  failure  time  axis  by  an  amount  cK  ;  then  the  adjusted 
points  can  be  connected  by  a  straight  line. 

Estimates  of  the  parameters^  and  are  obtained  directly  from  the 
straight  line  plot.  As  will  be  shown  below,  the  slope  of  the  straight 
line  is  ,  and  the  intercept  of  the  straight  line  on  the  vertical  axis 
is  -y^ln-0-,  from  which  one  can  obtain  the  value  of  *  In  general 
practice,  however,  it  is  more  common  to  obtain  the  value  of  ^-merely  by 
evaluating 


=  1 

where  t  is  the  value  of  failure  time  on  the  Weibull  plot  corresponding 
to  F(t)  =  0.632.  The  value  0.632  is  derived  from  the  fact  that  for  any 
value  of  0  ,  the  value  of  F(t)  obtained  from  equation  (l)  when  t/&-  =  1  is 
always  O.632.  In  this  study,  the  latter  method  for  determining  was  not 
used.  It  was  found  to  be  equally  expeditious  in  the  computer  program  to 
evaluate -Qr  by  equating  the  intercept  on  the  vertical  axis  to  -  In  -£h 
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To  show  that  the  slope  of  the  straight  line  plot  on  Weibull  probability- 
paper  is  ft  ,  we  start  with  equation  (l): 


F(t) 


and  obtain  the  following: 


1  _  tvf 

i  -  F(t) 


in  In  ^ 

1  )  =  i’lnt  -  ft  In  ■& 

rrmj  r  F 

(a) 

Equation  (ii)  is  seen  to  be 

of  the  form  Y  -  X  +  B 

when  we  let  Y  =  In  In  1 

1  ) 

(5) 

u 

X  -  In  t 

(6) 

B  *  -ft  In  ■&, 

(7) 

in  which  case  ft  is  the  slope  of  a  straight  line  plot  on  a  coordinate  system 
where  the  horizontal  axis  is  graduated  according  to  X  =  In  t  and  the 
vertical  axis  is  graduated  according  to  Y  =  In  In  (l/(l  -  F(t))).  Weibull 
probability  paper  has  its  axes  so  graduated,  consequently^  is  the  slope 
of  a  straight  line  plot  on  Weibull  probability  paper  of  F(t)  vs.  failure 
time,  t.  Note  that  the  intercept  of  the  straight  line  plot  is  B  and  that 
it  equals  -ft  In-#-.  It  is  seen  then  that-0-  can  be  calculated  once  values 
for  ft  and  B  are  obtained. 


It  should  be  noted  that  because  estimators  such  as  mean  ranks  or  median 
ranks  are  used  in  this  method  as  estimations  for  F(t),  the  resulting 
calculated  values  of ft  and  -9~ are  estimates  of  ft  and  -9-  and  should  be 
labeled^  and  £  to  indicate  that  they  are  estimated  values. 

Sometimes  in  life  testing,  the  test  is  cut  short  or  stopped  prior  to 
having  all  of  the  test  samples  fail.  For  example,  if  hO%  of  the  test 
samples  have  not  yet  failed  when  the  test  is  terminated,  the  test  is 
termed  a  1x0%  suspended  data  test.  If  60%  of  the  test  samples  have  not 
failed  by  the  time  the  test  is  terminated,  the  test  is  termed  a  60% 
suspended  data  test.  Suspended  data  tests  are  intuitively  expected  to 
be  less  accurate  than  non-suspended  data  tests  for  estimating  the  Weibull 
parameters  ft  and-#-  because  they  result  in  fewer  points  on  the  Weibull 
probability  paper  through  which  to  draw  the  best  fit  straight  line. 
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APPENDIX  II 


COMPUTER  PROGRAM 


W  £  I B  U  L V  C  M  6  0  0  0  0  »  T  5  0  0  .  . 

COMMENT. (TAA-102* 000000) ^CATALANO 

FTN . '  '  '  . 

LGO. 

0000000 0000 000 0000 00 00  '  -  —  ~ 

PROGRAM  PR685S ( INPUT  *  OUTPUT) 

C  '  '  " . 

C  SIMULATION  OF  WEIBULL  FAILURES 

C . COMPARISON  "OF'  MEDIAN  RANKS  AND  MEAN  RANKS 

C  AS  WEIBULL  RANKING  APPROXIMATORS 

C .  COMPARISON  OF  APPROXIMATIONS  WHEN  USING  "  40* » 60%’  ’AND  09TSUSPENSION 

C 

C  .  ***MA IN  PROGRAM***'  .  '  '  . 

C 

D I  MENS" I  ON  V (1 6  0  )  *  T  ( 160  )  v  RMEO  <  1 60  )♦  RMEN  ( 160  )  »  . 

SBETMED ( 20 ♦ 3 ) tFTMED <20 *3) * THMED (20 • 3> » 

. SbETMEN ( 20  «  3 ) » FTMEN ( 20 ♦ 3  ) . THMEN ( 20 » 3 ) V  .  . 

SX ( 1 60 ) *XRMEO(160) ♦ XRMEN (160) 

.  REAL  MDBMEO  (3)  ♦MD3MEN  (3)  »MDTMED  (3) VMDTMEN  <  3) . "  . .  '' 

C 

~C  '  PRTnT  HEADINGS 

c 

- PRINT  1 4  . 

J.4  FORMAT  (^1F»1X.17H**»*VARIABLES****»17X»8SH***********RESULTS - ST  A 

“SNDARD  ERROR  (FOR  20  MONTE-CARLO' SIMULATION  TRIALS)*** *******»//» 
S38X»24H******MEDIAN  RANKS****** » 32X » 24H*******ME AN  RANKS*******  » 
^PERCENT  A GE~ SUSPENSION  t  3SX  ♦  21HPERCENTA6E  "SUSPEnSI  ONV 
__S//*_1X*6HSAMPLE*22X*3H60%»16X*3H40%«  15X»4H  0%  »15Xf3H60‘i>*l6X» 

S3H40 % ♦  1 5X  » 4H  0%  »/»lXt4HSIZE»2X94HBETA»3X*5HTHETA»5X*4H8ETA»5X* 
S5HTHETA,SX,4H8ETA,5X»SHTHETA»5Xv4HBLTA»3X»SHTHETA,  fX ♦ 4HBE T A  * SX * 

”  SSHTHET A » SX *4HBET A»5xV5hTHETA»SX«4H8E7A*3X* 5HT HET A ) 

c 

C  '  "ASSIGNMENT  'OF  PARAMETER  "VA1ZOES" '(THCTA*“ BETA,  N) 

C 

DO"Y2~  I'2'=l  *6“ - - —  . . . . - . 

PRINT20 

. . "20  FORMAT!/)  '  ~  .  .  . ““  . 

THETA=5Q00.0*2**(I2-1) 

D0  -10  L'=i  »5 - ~ - — - — - 

PRINT  18 

"18"  FORMA  re/)  '  —  — - - 

BETA=(2**L)/4.0 

"""" .  DO  S  K=  1 » 6  .  ~ . 

N=S*2** (K-l ) 

'DO  6  M=l*20  -----  -  -  .  -  . — 
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o  o  o;o  o|o  !  n  o  r>  o  r>  j  !  ooo  ooo 


GENERATE"  A  SET  OF~N" RANDOM  ^NUMBERS 


DO  31  =  1  *N 
30  Y  ( I ) =RANF ( 0 ) 

IF  ( Y  ( I )  .LT.0.0)  GOOff 
3  CONTINUE 


RANK  ORDER  THE  SET  OF  RANDOM  NUMBERS 


DO  5  1  =  1  *  N 
1)0  7  J=I,N 
IFCY(I)-Y(J) >7*7*9 

9  S=Y  ( I)  “  . . 

Y ( I ) =Y ( J> 

Y ( J)=S 
7  CONTINUE 
5  CONTINUE 


GENERATE  N  MONTE-CARLO  FAILURE  TIMES 
N  MEDIAN  RANKS 
N  MEAN  RANKS 


. DO  11  1  =  1 *N 

T(I)  =  (THETA)*(-ALOG(1.'0-Y(I>  >  >**(1.0/BETA> 

KMED ( I ) = ( I -0 • 3) / (N+0 .4) . 

RMENd  )=I/(N*1.0) _ _ _ _ _____ 

11  CONTINUE 

-  -~N  ££  —  g— ^TRANSFORMED"  “A  XES  VMTuES  OF 

FAILURE  TIMES :X 

. MEDIAN  RANKS':  “XRMED 

MEAN  RANKS:  XRMEN 


DO  1  1  =  1 »N  _  ' _ 

X  d)=ALOG(T ( I )  ) 

XRMED ( I ) =ALOG ( ALOG ( 1 . 0/ ( 1 . 0-RMED ( I ) > ) ) 
1  XRMEN(I)=ALOG  (ALOG  (1.0/  d.O-RMENdf)  )  ) 
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nrin  j  |  !  o;nn{  |  j  jO  n  n 


perform  least  squares  fitting  for  median  ranks 


DO  40  1=1*3 
N=5*2** (K-l ) 

N=  I F  I X  (  (0.2*0 .2  *  2  **  (  1-1 7  T*'NT  ' 

Call  lsfit (n»x*xrmeq9A»b»c*d»e*F) 

BETMED (M* I )= (C«E-B*F)/(A*E-H*D>  _ 

FTMEDtM*  I)=(A*F-C#D>/(A*E-B*D) 

ThMED ( M » I  ) =EXP (- (FTMED ( M  *  I > ) /  ( BET MED (MVf)  )  ) 


PERFORM-LEAST-SQUARES-FITTING-FOR-MEAN-RANKS 


CALL  LSFIT (N »X * XRMEN* A*B*C*D*E *F) 

BE  T  M  E  N  ( M  i T>  ="(  C  *  E  -  B '* F )  /  <  A  * E  - d  *  D 
FTMEN  ( M »  I)  =  (A*F-C*D)/(A<*E-BttD) 

ThMEN(M» I)=EXP(-<FTMEN(M»I) )/(BETMEN(MV1) > > 
40  CONTINUE 

6  CONTINUE  . . 


. . .  '  PERFORM  'STANDARD'  ERROR  CALCULATIONS  .  ‘  '  " 

DO  50  1  =  1  »3  “ 

CALL  MDEV(I*BETA*BETMED*MDBMED) 

"  “  CALL  MDEV ( 1 1  BETA  * BETMENTMOBMEN) 

CALL  MDEV ( I ♦ THET A ♦ THMED ♦ MOTMEO ) 

XA  L  L ' M  D  E  V ( I ♦ T  H  £  T  A VT  H M  E  N » M  D  T  M  E N )  '  • 

50  CONTINUE 

'  print i 6»  n  *beta  »  thetavmdbmed  n evmdt med  m~»  mdbmedt2)~ ♦mdtmed ( 2 )  *mdbme 

SD(3) ♦MDTMED(3) *  MDBMEN ( 1 ) ♦MQTMEN(l) *  MDBMEN ( 2 ) ♦ MDTMEN ( 2 ) ♦ MDBMEN ( 3 ) » 

SMDTMEN  13)  -  - - - - - 

16  FORMAT (lX»I3»3X*F3.1*2X*F8.1t3X*F6.3*lX»FlQ.l»2X*F6.3«lX*FI0.1 » 

S~2  X' » F5  •  3  *  1  XTF"9TTV3X*"F  6"»"3  *1 X ,  FT07IT2xVF6.3, 1X»F  1 07TV2X» 75 . 3 ♦ 2X,F9 .1) 
8  CONTINUE 

TO  CONTINUE  . . 

12  CONTINUE 
STOP 
END 


3k 


o  no  j  i  |  o  o  o  o  r>  o 


SUBROUT  INES***' 


SUBROUTINE  FOR  LEAST  SQUARES  FITTING 

SUBROUTINE  LSFIT (N*T»R*A*B*C*D*E»F> 
kEAL  T ( 160 )  * R ( 160 ) 

SUMX  =  0  •  0 
DO  1  J= 1  *  N 
1  SUMX=SUMX*T ( J) 

"  A * SUM X 

d=N/1.0 
SUM  Y  =  0""*0 
UO  2  J=1*N 
2SUMY=SUMY+R(J) 

C=SUM Y 

. S  U  M  X  2  =  0  •  0 

UO  3  J=1 *N 

3  SUMX2=SUMX2+T ( J) **2 
U=SUMX2 
E  =  A 

SUMX Y=0 • 0 
UO  4  J=1,N 

4  SUMXY  =  SUMXY *T (J)*R(J) 

(■  =SUMXY  "  .  . 

RETURN 
end 

s u b r outTne'for^lculatTng  s ta'noari (Terror 

SUBROUT INE  MOEV ( I »KNO*N»EST ♦MENDEvT 

REAL  KNOWN»MENDE_V  (3)  »EST (20  >3)  >DEV2  (20*3) _ 

UO  1  M=i *20 

1  DEV2 ( M » I )  =  ( KNOWN-EST  ( M  *_I_)  )  »»2 _ 

SuHQEV  =  0 . 0  '  ’  "  . .  ~ . . 

DO  2  M= I » 20 

2  SUNDEV=SUMDEVfDEV2(M* I) 

UUMMY=SUMDEV/20.0 
MENDEV  ( I )  =SQRT  (OTJMMYT 
RETURN 

END . .  '  ‘ 

0 00 0 00 U 00 00 0000 0000000 
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APPENDIX  III 

CHARTED  PRINTOUT  OF  RESULTS 
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CHART  NO.  1 

PRINTOUT  OF  STANDARD  1RROR  WHEN  USING  MEAN  RANKS 
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STANDARD  ERROR 


RANKS******* 


*e#-u-VARIA6LES**** 

PERCENTAGE 

SUSPENSION 

sample 

60  it 

40% 

0% 

SIZE"  ‘ 

bETA 

THETA 

r*": 

THETA 

BETA 

THETA 

BETA 

TFHET  A 

5 

.5 

5000.0 

*7.469 

4793735.8 

*7.469 

4793735.8 

.234 

8128.7 

10 

,b 

5000.0 

. .305 

195047.7 

.232 

46925.9 

.130 

4495.6 

20 

•  b 

5000.0 

.241 

58775.4 

.175 

14709.9 

.111 

2936.5 

40 

.  5 

5000.0 

- .208 

95934.9 

.150 

18468.0 

.008 

2  1 94.6 

80 

•  s 

5000.0 

.136 

1 3663 . 8 

.112 

4114,4 

.  059 

9  6  4  . 

160 

c  b 

5000.0 

‘  .089 

6871.7 

.070 

3040.2 

.04  6 

9  6  8  .  e. 

5 

1  .0 

5000.0 

3.974 

28086.3 

3.974 

28086.3 

.8  77 

2295.1 

....  1Q  - 

1.0  “ 

'5000.  0| 

.619 

.  8929.7 

.483 

6144.2 

.285 

1804  .  v 

?0 

1  .0 

5000.01 

.36  8 

13606.8 

.307 

4681.5 

.  2 1  0 

1  0  2  2  .  : 

4  0 

i  .0 

50  00 . 0 | 

.317 

2433.2 

.223 

1 435 . 3 

.163 

81  7.3 

eo 

1 .0 

50  0  0 . 0  j 

.186 

2890.1 

.143 

1370.8 

.089 

57  7  • 

160 

1.0 

b  0  0  0  •  0  I 

i 

.158 

254  1,2 

.134 

1139.6 

.092 

4  36 . 7 

5 

2.0 

5000.0 

5.793 

11001.4 

5.793 

11001.4 

1.002 

1385.1 

10 

2 . 0 

5000.0 

2.967 

6762.9 

1.220 

4010.4 

.513 

7  58.3 

20 

8.0 

5000.0 

!  .814 

3033.3 

•  66  1 

1637.0 

.460 

589.  l 

40 

2, 0 

5000.0 

;  .609 

2160.4 

.555 

1164.1 

.389 

474.2 

80 

2.0 

5000. 0 

.369 

954.3 

.335 

602.0 

.240 

247.3 

'  160 

8 . 0 

. 50  0  0.0 

! . ".4  1 6 

1 

1156.9 

.344 

536.3 

-.24  0 

192.5 

5 

4 , 0 

5000.0 

21  .403 

1701.4 

21 .403 

1701.4 

2.766 

591 .3 

10 

4.0 

"  5000.0 

!  ■'  2.641 

2330.0 

2.443 

1593.2 

1.6  70 

414.7 

20 

4.0 

5000.0 

1.890 

1661 .8 

1  .467 

929.9 

1.003 

272.1 

40 

4.0 

. 5000.0 

. -  1 .20  3 

. 607.1 

. -.983 

- 348.2 

-.797 

. -  20  5.9 

80 

4.0 

5000.0 

|  1.064 

623.6 

.952 

360.4 

.594 

1  5  8  , 

160 

4 . 0 

5000.0 

-  '  .624 

.  300.8 

.513 

187.8 

.318 

-  93.3 

S 

6 . 0 

5000.0 

'  24.652 

811  .8 

24.652 

811.8 

2.852 

325.5 

10 

8.0  ' 

' 5000.0 

|  -  5.030 

. 1117.9 

4.452 

-  -  748.6 

2.559 

-  2  5  f~.  #  j-; 

20 

.  0 

5000.0 

1  2.920 

650. 1 

2.356 

421.2 

1.728 

206.9 

40 

r‘ .  0 

5000.0 

i  2.678 

425.4 

2.556 

267.9 

1.497 

121.0 

80 

6 . 0 

5000.0 

!  1.615 

291.8 

1  .516 

165.5 

1.0)6 

61.4 

160 

-.0 

5000.0 

j  1 .459 

234.3 

1.101 

133.7 

.  b4  7 

53.2 
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STANDARD  ERROR 


f  •  . .  ‘ 

#####*#MEAN  RANKS******* 


****VARI  A6LES**** 

PERCENTAGE  SUSPENSION 

SAMPLE 

6011 

40% 

0% 

SIZE 

BETA 

THETA  j 

BETA 

THETA 

BETA 

THETA 

BETA 

THET  A 

5 

.5 

10000.0 

2.629 

#4440491 . 8 

2.629 

*4440491 .8 

.354 

21169.2 

10 

.5 

10000.0 

.819 

*0765708.3 

.458 

22425170.1 

-  .195 

12681  .u 

20 

.5 

10000.0 

.244 

460722.5 

.179 

85459.1 

.116 

8162.6 

-  40 

— T5T 

10000 .0 

•  1 46 

37054.1 

“  .119 

12070.1 

.0  72 

3669.2 

60 

.5 

10000.0 

.120 

12775.7 

.095 

7090.5 

.  056 

2844.1 

160 

.5 

10  000. 0 

.075 

9613.8 

.061 

3814.8 

.044 

1741.9 

5 

1.0 

10000.0 

4.778 

38094.2 

4.778 

38094.2 

.  38  1 

5331.6 

'  TO . 

"t:o: . 

10000.0 

.861 

28821.0 

-  .68  7 

1 6624 . 8 

302 

4504.4 

‘  20  « 

1.0 

10000.0 

.362 

4  7  6  6  3 . 2 

.309 

16022.6 

.202 

284 6.4 

4  0 

”1.0 

10000.0 

.310 

17122.1 

.277 

7238.7 

-  .190 

2806.7 

80 

1 .0 

10000.0 

.187 

5040.9 

.157 

2677.9 

.103 

llbO.v 

1 60 

1.0 

1000010 

.173 

4572.9 

.133 

2402.2 

.0  75 

926.4 

5 

2.0 

10000.0 

3.183 

9223.4 

3.183 

6223.4 

.982 

2193.7 

‘  10 

"2 . 0 . 

”1  0000.0 

.  8  7  6 

8515.8 

.826 

—  ■  5293.5 

92 

1597.1 

20 

2.0 

10000.0 

.613 

4009.6 

.572 

2859.0 

.43  7 

1241.3 

4  0  . 

”210 

10000.0 

.506 

2498.5 

-  .362 

1292.1 

.269 

625 . 3 

80 

2.0 

10000.0 

.361 

3702.0 

.287 

1814.6 

.205 

797.9 

160 

2.0 

10000.0 

.314 

2011.4 

.251 

1101.7 

.158 

511.5 

5 

4.0 

1  0000. 0 

17.023 

3178.2 

17.023 

3178.2 

3.211 

1334.5 

10 

4.0 

10000.0 

2.460 

4279.8 

2.176 

2915.6 

1.470 

844  .  i 

20 

4.0 

10000.0 

2.653 

3229.4 

1.805 

1918.4 

1.016 

627 . 1 

40 

** .  0 

10000.0 

1.250 

2449.4 

1.108 

1430.1 

.783 

438 . 6 

60 

4.0 

10000.0 

.611 

674.2 

.508 

378.3 

.362 

262. U 

160 

H.o 

10000.0 

.620 

666.2 

.534 

447.6 

.316 

226.6 

5 

8.0 

10000.0 

19.457 

1994.3 

19.457 

1994.3 

2.932 

625 . 7 

10 

8.0 

10000.0 

3.908 

2276.4 

3.366 

1568.3 

2 . 1  D  ( 

462.4 

20 

8 . 0 

10000.0 

3.690 

1576.6 

3.036 

916.2 

2.137. 

250.6 

40 

8.0 

10000.0 

2.437 

917.3 

2.248 

548.2 

1.659 

214.5 

80 

6.0 

10000.0 

1.501 

528.8 

1.193 

291 . 1 

•  7d8 

132.0 

160 

8.0 

10000.0 

1.056 

305.1 

.879 

197.1 

"  .603  ' 

109.5 

39 
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« ■&#•&« * ■» M E A N  RANKS******* 


****  VAR  TABLES***'* 

PERCENTAGE 

SUSPENSION 

SAMPLE 

60S 

40% 

0% 

SIZE 

bET  A 

THETA 

BETA 

THETA 

BETA 

THETA 

BETA" 

'  'THETA 

5 

.5 

•20000.0 

4.288 

*3955750.0 

4.288  *3955750.0 

.272 

34192.8 

10 

.5 

20000.0 

.323 

150598 18.5 

.325 

1507381.4 

.134 

13472.2 

20 

.5 

20000.0 

.207 

1186650.0 

.171 

139793.8 

.127 

1  1376.  v 

40 

•  i5 

20000.0 

.141 

62708. 1 

"  '  .  11  4 

32487. 1  ‘ 

"  .075 

1  1377.0 

80 

.5 

20000.0 

.134 

1 16422.5 

.109 

25640.8 

.073 

5778. 6 

160 

.5 

20000.0 

.097 

31673.3 

.074 

11886.5 

.04  8 

3371.6 

5 

1.0 

20000.0 

7.635 

35196.9 

7.635 

35196.9 

.429 

9123.4 

10  , 

1.0 

20000.0 

"  .765 

. 99045.9 

.729 

50296.8 

.327 

'  10949.3 

20 

1  .0 

20000.0 

.357 

26047.9 

.34  0 

13189.7 

.241 

6365." 

40 

1 .0 

20000.0 

.302 

28543.5 

.257 

12235.7 

.187 

4428.6 

80 

1.0 

20000.0 

.190 

12713.7 

.162 

5945.1 

.122 

2519.6 

160 

1.0 

20000.0 

.150 

5980.9 

.117 

3398.3 

.077 

1525.6 

5 

2.0 

20000.0 

6.350 

72237.7 

6.350 

72237.7 

1.092 

6098. 1 

'  io'" 

‘  2  <  0 

20000.0 

1.259 

16059.0 

1.188 

11344.2 

•  666 

'  3627 . 

20 

2.0 

20000.0 

.690 

11966.2 

.604 

7079.0 

.453 

2438.8 

""40 

2.0“ 

■'"20  0  0  0.0 

.536 

7279.0 

.479  . 

3365.6 

.331 

1632.3 

80 

2.0 

20000.0 

.  3b0 

4132.0 

.311 

2484.7 

.219 

1254.2 

160 

2.0 

20000.0 

.260 

2954.9 

.199  ' 

1907.7  ~ 

•  1^7 

"1006.6 

5 

8.0 

20000.0 

’  8 . 753 

8710.7 

8.753 

8710.7 

1  .44  1 

2431.1 

10 

4.0 

20000.0 

1.984 

13555.8 

1  .857 

8449.5 

1.413 

1952.  v 

20 

44.0 

20000.0 

2.124 

5204.6 

1.698 

'  3100.5 

.978 

■1368.* 

40 

44.0 

20000.0 

1.233 

6249.9 

1.042 

3053.5 

.717 

809. <4 

80 

4 . 0 

20000.0 

;  . . .967 

. 3149.9 

'  .811 

"“1 705.8 

.583 

-  519.2 

160 

44.0 

20000.0 

|  .604 

1713.6 

.439 

964.1 

.266 

465  •  7 

~ . 5" 

8  ~0 

"  '2000  0VO 

1 

1 

\  *3.228 

4002.7 

*3.228 

-  4002.7 

4.04  7 

1367.7 

10 

8.0 

20000.0 

!  5.115 

4376.4 

4.205 

3168.6 

2.463 

88  0.2 

“  20 

6,0 

20000.0 

3.711 

3430.0 

3.38V 

2022.6 

2.147 

64  0  .  h 

40 

6 . 0 

20000.0 

2.181 

1389.6 

1.881 

865.8 

1  .457 

343.  - 

80 

6.0 

"20000.0 

1.981 

1192.8 

1  .522 

671.3 

1,10  7 

3  72.8 

160 

8.0 

20000.0 

1.209 

681.0 

1.011 

340.0 

.686 

1  4  8 . 1 

1|0 

-i  rv 


STANDARD  ERROR 


f  '  • .  . .  '  . .  '  '  . .  ’ 

-fr-a-ME AN  RANKS******* 


***•«•  variables**** 

PERCENTAGE 

SUSPENSION 

sample 

60S 

40% 

0% 

size 

bETA 

THETA 

BETA 

THETA 

BETA 

THETA 

BETA 

"THETA 

t 

f 

5 

.5 

40000.0 

!  . -.520 

1  7076942.5 

.520 

1  7 0 7 6942.6  ‘ 

193 

61540.3 

10 

.5 

40000.0 

.371 

2951130.2 

.315 

736304.9 

.156 

39961.2 

20 

.5 

’40000.0 

.30? 

849440.1 

.239 

1  84436.8 

.092 

25933.5 

40 

.5 

40000.0 

.17? 

838168.0 

.153 

1 32430 • 2 

.096 

11940.2 

80 

.  5 

40000.0 

- .095 

-  39103.1 

.073 

2017  7.7 

.049 

1  0402. 1 

160 

.5 

40000 . 0 

.076 

43691.2 

.059 

■- 

21528.9 

.042 

7918.1 

5 

1.0  ’ 

400  00 . 0 

2.765 

174940.0 

2.765 

174940.0 

.4  09 

20646.0 

10 

1  .0 

40000.0 

.628 

106860.0 

.525 

43679.2 

.310 

9687.4 

. . 20 

l.o . 

40000.0 

.417 

93283.4 

. .331 

25918.8 

.200 

9648.7 

40 

,  1.0 

40000.0 

.287 

151932.9 

.245 

38834.8 

.180 

920 8 . 1 

80 

1.0 

40000.0 

.207 

19031.1 

.169 

8365.9 

.125 

3929.9 

160 

1.0 

40000.0  : 

.167 

13987.3 

.137 

7529.8 

.086 

3169.2 

-s_ 

2".  O' . 

40000 . 0 

i 

1 

9.450 

35979.0 

9.450 

35979.0 

1 .369 

-  10327.1 

10 

2.0 

40000.0 

1.119 

29976.4 

1.084 

19052.3 

.705 

6845 . 9 

20 

"’2.0 

'40000.0 

1.089 

19473.4 

'  .623 

10778.1 

.3  7b 

4506.2 

40 

2.0 

40000 .0 

.633 

9254.0 

.588 

5451.3 

.328 

24l6 . 2 

. 80 

"’2.0 

”40  0  0  0.0 

.437 

12158.4 

.375 

5487.9 

.270 

2568 . 2 

160 

2.0 

40000.0 

.349 

5880.7 

.299 

3525.3 

.201 

1642.8 

5- 

4'ro' 

t 

”40000.0 

7.333 

.  57435.7 

7.333 

57435.7 

1.92b 

4740.6 

10 

4.0 

40000.0 

1.766 

20490.2 

1.594 

1 3809.9 

1.161 

3994.3 

"  20 

”470'  ' 

"4  0  0  0  0.0 

1.518 

9740.0 

1.290 

5244.7 

.901 

1  56b . 9 

40 

4.0 

40000 . 0 

1.382 

4956.6 

1.254 

3019.7 

.  b93 

1224.3 

80’ 

4  70 

40000.0 

- .80  7 

- 4488,3 

.664 

-  2715.1  - 

-.45  0 

1176. 3 

160 

4.0 

40000 .0 

.596 

3790.5 

.470 

2066.2 

.306 

f  5 1  •  6 

5 

8  ”  0” 

40000.0 

10 

8.0 

40000 . 0 

”20” 

8.0” 

— 40000.0 

40 

8.0 

40000.0 

80 

8  30 

40000.0 

160 

5.0 

40000.0 

18.391  - 4704.8 

6 . 753  5783.2 

3.348  - 4030.3 

1.863  2334.3 

1.317 - 1773.6 

1.202  1296.8 


18.391  4704.8 

4.497  4186.8 

2.525  -  2710.-3 

1.720  1547.1 

1.012  . ~  1  136.9 

.906  745.3 


4.655 
1.963  1743.  o 

~ 1 . 64  1  .  1254 . 1 

1.370  925.4 

.683  “ 57  3.5 

•608  335 ; a 


c>  a- 


STANDARD  ERROR 


*#**«*» ME AN  RANKS******* 


e»««yARI  A6LES***'* 

PERCENTAGE  SUSPENSION 

sample 

SOS  40%  0% 

SIZE  bETA  THETA 

BETA  THETA  BETA  ThETA  BETA  THETA 

l! 


5~ 

--g— 

80000.0 ! 

5 . 854 

*5670590.2 

5.854 

*5670590.2 

.215 

105767.3 

10 

.5 

80000.0 ; 

.215 

*  4  3  4 7  3  6  4 .6 

.189 

1098^985.9 

.  1  4b 

5  t>  5  6  2  •  _• 

“2D 

- V5~ 

. 80000 . 0 | 

.197 

1865397.8 

.  1  64 

392070.4 

.104 

66366.6 

40 

.5 

80000 .0  j 

.137 

162364.9 

.119 

66025.1 

.0  77 

2729b .  9 

80“ 

- -5- 

800  00.0  I 

-  .103 

122425.3 

.088 

56 116.0 

.  Ob8 

24966 . 5 

160 

.5 

80000.0 ' 

.....  1 

.060 

46840.0 

.04  6 

25320.1 

.031 

1 2566 . 6 

•  5’ 

■"  1 .  O'" 

80000.0 

.632 

374613.7 

.632 

374813.7 

.7)0 

3/646.  b 

10 

1.0 

80000.0 

•  508 

876 3 69. 9 

.  50a 

328428.0 

.331 

3073‘+ . 

"  20 

-  1 .0“ 

80 o o o ; o 

. .344 

152111 .8 

.286 

72163.8 

.226 

25326.8 

40 

*  :  1.0 

800  00 . 0 

.442 

124317.4 

•  3  5  o 

4 6 6 04.2 

-  i  v6 

16157.. 

80 

. 1.0 

'  80000.0 

.250 

53494.5 

.199 

27259.8 

.  i25 

1  3  6  o  6  •  / 

160 

1.0 

80000 . 0 

.206 

69729.3 

.168 

?6738 . 4 

.112 

5966 .  / 

5 

2.0 

80000.0 

22*993 

242446.1 

22.993 

24 V4 98. 1 

.79  1 

19650.7 

10 

2.0 

80000.0 

1.063 

75930 . 8 

.839 

45166.3 

.4/0 

14370 . 1 

20 

2.0' 

80  0  0  0 .0 

.793 

423 26.8 

.625 

20258.4 

.  444 

7861.6 

40 

2.0 

80000.0 

.713 

37945.0 

.615 

19253.0 

.401 

6  6  3  9  .  t 

80 

2.0 

"  80  00  0 . 0 

.5  1 6 

27514.6 

.427 

14068.9 

•  2bO 

4622.7 

160 

2 . 0 

8  0  0  00.0 

.416 

16384.4 

.350 

9845.7 

.219 

2549.4 

. . 5 

4.0 

*60000.0 

39.543 

“50023.1 

39.543 

50023.1 

1.697 

9544.6 

10- 

4 . 0 

80000.0 

1.721 

46515.8 

1.635 

2  Q  2  2  6 . 7 

1.154 

7097.2 

20 

4.0 

80000.0 

1 .349 

13414.4 

1.081 

8519.5 

.696 

3379.6 

40 

4 . 0 

80000.0 

1.017 

1 7277.6 

.631 

8843.8 

.659 

3106.  i 

80 

4.0 

*60  0  0  0.0 

.863 

8733.4 

.745 

5178.6 

.894 

2671 .0 

160 

l 

1 

!* 

.  • 

80000.0 

.606 

4942.7 

.464 

2712.7 

.305 

1705.4 

5 

e .  0 

80000.0 

'  24.269 

15204.1 

24.269 

15204.1 

3.432 

6457 .6 

1  0 

4  .  o 

80000.0 

5.74? 

10399.2 

5.659 

6759.1 

2.654 

2944 .a 

--  2Q- 

h.  O’ 

80  0  0  0. 0 

3.025 

10930.6 

2.426 

6537.7 

"  1.593 

"  2107.5 

40 

8.0 

80000.0 

1.829 

4131.7 

1  .630 

2742.0 

1  .  1  -»5 

1486.0 

RO 

H  _  fl 

8  0  0  0  0 . 0 

1.188 

4110.8 

1.009 

2569.9 

.  645 

1252.1 

160 

8.0 

80000 .0 

1.572 

3383.8 

1.238 

1512.4 

.608 

690 . 1 

STANDARD  ERROR 


( 

RANKS*** **** 


****VARIABLES**** 

PERCENTAGE  SUSPENSION 

SAMP 

SIZE 

LE 

60S 

40% 

0% 

BETA 

THETA 

BETA 

THETA 

Beta 

theta 

'BETA 

THETA 

5 

.5 

160000.0 

! 

2.601 

*3324794.5 

2.601 

*3324794.5 

.167 

235484.4 

10 

.5 

160000.0 

.402 

5504098 • 1 

-  .277 

1521504.1 

.191 

102813.7 

20 

.5 

160000.0 

.241 

1 1543474.9 

.195 

1 36^642 . 8 

.153 

113046.8 

^0 

.5 

160000.0 

•  1  86 

1943901.0 

.148 

433706.3 

.093 

80017.4 

80 

•  s 

160000.0 

.128 

”769383.3 

.086 

1 6665 1.6 

.049 

43660.4 

160 

.  5 

160000.0 

.062 

86914.8 

.055 

44124.8 

.044 

27758.4 

" . 5 

1.0 

160000.0 

"  7.489 

542124.1 

7.489 

548124.1 

•  56  3 

84344.2 

10 

1.0 

160000.0 

2.018 

237742.5 

1.127 

131731  .7 

.317 

517S6.1 

20  v 

«  l.o" 

160000.0 

'  ".538 

547686.4 

”  .352 

171 842.4 

- .236 

"—52  792.8 

4  0 

1.0 

160 00 0.0 

.263 

1 24685 . 3 

.196 

55778 . 4 

.129 

17962. 0 

80 

T  .  0 

160 000.0 

.225 

86972.5 

.157 

46806.0 

'  -  -■  1 1 4 

22425.5 

160 

1.0 

160000.0 

.136 

57099.6 

.102 

28077.8 

.  066 

13318.2 

‘”"5“ 

16000  0.0 

54.392 

- 76935.5 

54.392 

.  76935.5 

.708 

35594.5 

10 

2.0 

160000.0 

1.198 

132753.7 

1.108 

77820.2 

.580 

25056 • 3 

20 

. 2.0" 

TbOOOO.O 

- . 1.014 

205919.7 

.879 

52837.2 

.461 

"24294.3 

40 

2.0 

160000.0 

.607 

65914.0 

.528 

37286.9 

.359 

11505.2 

. 80 

2.0  ' 

1 60  0  0 0 . 0 

- .391 

”33629.2 

.319 

18695.5 

.213 

- 7794.3 

_  160 

2.0 

160000.0 

.339 

31807.9 

.274 

17383.0 

.  1  7  4 

662 1.3 

5- 

4.0 

160000.0 

"22.695 

~ 44114.6 

22.695 

44114.6 

1.948 

21968.5 

10 

4.0 

160000.0 

2.200 

43813.5 

2.126 

33752.8 

.915 

1 1574. b 

"20 

•"'4.0"" 

160000.0 

. 2.371 

31267.7 

1  .800 

20013.9 

1.017 

9651.5 

40 

4.0 

160000.0 

.992 

28297.3 

.840 

16263.3 

.600 

8020.1 

4.  O' 

160000 .0 

—  I .001 

- 19599.4 

■  '•  ;  740 

- q278#1 

.45 1 

- 4480.6 

160 

4.0 

160000.0 

.557 

1 1405.3 

.394 

657  8.4 

.281 

3309.4 

. 44.571 

26029.1 

44.571 

5.050 

_  S 

35.0 

160000.0 

• 1 

O  l)  1  D  •  D 

10 

6 . 0 

160000.0 

7.859 

24018.4 

3.698 

17887.0 

2.430 

52  7 3 • 9 

20 

d.O" 

160000.0 

3".  499 

22257.6 

”3.050 

- 13341.0 

2.032 

4504.3 

40 

8.0 

lbOOOO.O 

2.654 

12774.6 

2.401 

7037.1 

1.481 

3168. b 

80 

'6.0 

160000.0 

1 .796 

9719.3 

1.453 

■  ~  5937.2 

‘  .867 

2915.7 

1  60 

6.0 

160000.0 

1 .450 

7063.0 

1.201 

4097.7 

.812 

1804.4 

CHART  NO.  ? 

PRINTOUT  OF  STANDARD  PRROR  WHHN  USING  MFDIAN  RANKS 


bh 


STANDARD  ERROR 


RANKS******* 


**#*VAR TABLES**** 

PERCENTAGE 

SUSPENSION 

sample 

SIZE  BET  A 

THETA  ' 

60% 

BETA  . THETA 

40% 

*  BETA  THETA 

0% 

— BETA  THETA 

( 

5 

•  s 

5000.0 

*9.794 

1217960.9 

*9 , 794 

1217960.9 

.2b4 

6623.4 

10 

.  5 

5000.0 

.381 

8/063.5 

.276 

27441.6 

.121 

4050.5 

20 

.5 

5000.0 

.261 

32312.8 

.175 

10580.0 

.100 

2711.5 

AO 

.5 

5000.0 

.  224 

59936.4 

.154  - 

1 4356 . 9 

■'  .  08  1  ~ 

2628.3 

80 

.5 

5000.0 

.142 

95 13.8 

.113 

3368.8 

.055 

95  0.3 

160 

.5 

5000.0 

.  0  8m- 

5561  .6 

.  066 

2  b  7  6.8 

.04  3 

945 . 2 

5 

1.0 

5000.0 

5 . 0  -8  d 

16566.4 

5.0  86 

1 6566 . 4 

1,054 

2142.6 

10 

1.0 

5000.0 

.766 

b3o  5.9 

.574  * 

4752 .7 

.295  * 

1711.6 

20 

1.0 

5000.0 

.  3oo 

9414.8 

.301 

360  / .  3 

.  1 1*  8 

1004,9 

40 

1.0 

5000.0 

.  3o2 

1925.3 

.244 

1252.6 

.168 

792.8 

80 

1.0 

5000.0 

.  1  8  0 

2335.1 

.135 

1200.4 

.0  79 

561.1 

160 

1.0 

5000.0 

.15  5 

2151.3 

.129 

1019.4* 

.086 

'  427.9 

5 

2.0 

5000.0 

7.447 

7650.9 

7.447 

7650.9 

1.177 

1312.1 

10 

2.0 

5000.0 

3.501 

5092.1 

1.396 

3210.3 

*  .502  “ 

"  696 . 3 

20 

2.0 

5000.0 

.87m 

2382.1 

.677 

1350.3 

8 

551.8 

40 

2.0 

"5000.0 

.609 

1746.1 

.564 

1025.6"  ' 

*  .371 

479.5 

80 

2.0 

5000.0 

•  33d 

816.8 

.314 

547.8 

.225 

248.1 

160 

2.0 

5000.0 

.379 

984.2 

.315 

468.2 

*  .223  ' 

192.9 

5 

4.0 

5000.0 

27.23b 

1463.8 

27.236 

1463.8 

3.383 

559.9 

TO 

4.0 

"5000.0 

3.137 

1883.2 

2.832  " 

1341.0 

1.764 

'  398.4 

20 

A  .  0 

5000.0 

2.06/ 

1372.0 

1.452 

805.0 

.853 

273.1 

40 

4.0 

5000.0 

1.201 

493.8 

.925 

301  .5 

.737 

205.1 

80 

4.0 

5000.0 

1.066 

526.8 

.960 

322.0 

.586 

15  7.0 

160 

4.0 

5000. 0 

.60  3 

263.2 

.492 

1/4.6 

.302 

.  94.8 

5 

8.0 

5000.0 

'  31.888 

697.4 

31.888 

697.4 

3.009 

328.8 

10  * 

r>  .  0 

50  0*0.0 

b .  9  6  2 

909.5 

5.032 

6<cT,  0 

2 .344 

. *  152.9 

20 

8.0 

5000.0 

2.677 

559.8 

2.220 

383.8 

1  .565 

208.2 

40“ 

—8.0 

— 50  0  0.0 

■  ~  2.654 

365.1 

2.528  - 

. . 238.3  - 

1 .364 

”118.2 

80 

8.0 

50  0  0.0 

1.620 

239.6 

1.375 

141.7 

.952 

60.3 

T60"~ 

— 6.D 

5000  •  0 

1  .437 

209.2 

1.061  - 

- 123.9  - 

”’.600 

- 52.2 

standard  error 


«««ae«ME f) JAN  RANKS****** 


#***VAR1 ABLES****  PERCENTAGE  SUSPENSION 


sample 

6054 

40$ 

0% 

SIZE 

BETA 

THETA 

BETA 

THETA 

'  BET  A'" 

THETA  — 

BETA’ 

‘  THET  A 

5 

.5 

10000.0 

3.344 

41443606. 3 

3.349 

41443606.3 

.488 

17636.0 

10 

.5 

10000.0 

.991 

<*7996729.5 

.535 

6707458.6 

.1  95 

10664.6 

20 

.5 

10000.0 

.2  74 

246833.4 

.  166 

59160,2 

.113 

7500.9 

40 

VS  ' 

10000.0 

.151 

24/10.7 

.125 

9380.7 

.  0  74 

3438 . 0 

80 

.5 

10000,0  ! 

.  1  34 

9750.2 

.  )  04 

6029,3 

.060 

2732.7 

160 

•  s 

looo o.o  ; 

.070 

7595.5 

.057 

3? 82.  1 

.04? 

1701.3 

5 

1.0 

1  0  0  0  0. 0 

6.07b 

20448.8 

6.076 

20448,8 

.  3  4  1 

4854. 7 

10 

'  1.0’" 

io ooo.o 

1 . 0  5  h 

19527.5 

.815 

12672.1 

'  .338 

4254.9 

20  " 

1.0 

1 0  0  n  0  .  u 

.  394 

320 1<  8.3 

.33  3 

124 73. 7 

.  2  0  3 

2  6  5  i  .  1 

40 

1.0 

10000.0 

.291 

12983.3 

.266 

6077.9 

.178 

2703 o 8 

80 

1*0 

10000.0  ! 

1  .174 

3988.7 

.148 

2292.5 

•  0  ->  6 

1123.5 

1 60  " 

1  .0 

10000.0  | 

j  .175 

3863.7 

.136 

2183. 7 

.077 

90  7.2 

5 

2.0 

10000.1) 

4.  124 

6953.5 

4.124 

6953.5 

1.10  7. 

2185.3 

10 

2 . 0 

1 0  0  00.0 

1.00  0 

6192.8 

.936 

4157.2 

.596 

1  5 1 2 . 6 

20 

2.0 

10000.0 

.606 

3211.3 

.532 

2431.5 

.38  0 

1180.0 

40 

2.0 

10000.0 

.813 

2009.9 

.353 

1117.6 

.25  7 

805.3 

80 

2.0 

10000.0 

.358 

3154.7 

.278 

1642.0 

.199 

783.0 

lbO 

2.0 

10000.0 

.321 

1  763.2 

.260 

1001.0 

.166 

494. 1 

5 

4.0 

10000.0 

j  21.750 

2925.3 

21.750 

2925.3 

3.653 

1346.8 

10 

4.0 

10000.0 

2 . 752 

331 1 .8 

2.330 

2392.0 

1.394 

860  •  5 

20' 

4 . 0 

10000.0 

j  3.30b 

2639.7 

2 . 0  0  9 

1654.4 

.974 

605.7 

40 

** .  0 

10  0  0  0. 0 

!  1.17c 

2089.7 

1.025 

1265.2 

.640 

415.9 

80 

4 , 0 

10000.0 

j  .557 

592.3 

•  454 

358.9 

.334 

265.4 

160 

4.0 

10000.0 

,  .614 

i 

610.4 

.  524 

425.1 

.30  0 

22  7.4 

5 

8.0 

10000.0 

r 

25.213 

1  763.6 

25.213 

1 763.6 

3.1  74 

652 . 4 

10 

8.0 

10000.0 

•  4.47  o 

1880.8 

3.635 

1331.3 

2.109 

470.8 

20 

>'.0 

loono.o 

!  3.839 

1293.7 

3.017 

771.7 

1.97  0 

239.0 

40 

8.0 

10000.0 

2.448 

80  0.3 

2.192 

4  95.0 

1  .Sbl 

214.3 

80 

8.0 

10000.0 

j  1.253 

476.4 

1.022 

250.8 

.6«2 

129.9 

160 

8.0 

10000.0 

;  .914 

251  .5 

.775 

173.7 

•  55  1 

107.6 

^6 


STANDARD  ERROR 


*»*«**MEDIAN  ranks****** 


****VARIABLES***» 

PERCENTAGE  SUSPENSION 

SAMPLE 

60S 

40% 

0% 

SIZE 

BETA 

THETA 

8ETA 

"THETA 

“BETA" 

THETA' 

BETA- 

THETA 

5 

.5 

20000.0 

5.409 

33645872.1 

5.409 

33645872, 1 

.30  1 

27622.1 

JO 

.5 

20000.0 

.391 

4430373.9 

.382 

645  7  0  8.1 

.117 

12516.6 

20 

.5 

20000.0 

.204 

41 2842.9 

.160 

79657.2 

.113 

1074^.5 

40 

•  B 

20000.0 

.  137 

57502.6 

.103 

26432.4 

.064 

10888.2 

80 

.5 

20  0  0  0. 0 

.  1  34 

606K 5 . 5 

,  1  06 

2090 1 .9 

.069 

5 4  6  8 . 0 

160 

.6 

2  0  0  0  0 . 0 

.  1  0  i 

25057 .8 

.  0  7  6 

10329.7 

.  04  V 

325o  .  4 

5 

1.0 

20000.0 

9.660 

22810.5 

9.660 

22810.5 

•  450 

8415.4 

. 10  .. 

1  .0' 

20000.0 

.920 

66594.0 

.  8So 

37966.6 

'  .332 

‘  10299.9 

20 

1.0 

20  0  0  0 . 0 

.3  74 

19506.2 

.362 

1  0  8 0  0.1 

.281 

6131.3 

40 

i.O 

20000.0 

.278 

21463.2 

.239 

10090.3 

.175 

7204.3 

80 

1.0 

200  00.0 

.  1  80 

10491 .9 

.156 

5260.4 

.117 

2479.6 

160 

1.0 

20000.0 

.  1  54 

5095.0 

.119 

' 3109.3 

.075 

1480.5 

5 

2.0 

20000.0 

8.154 

46489.2 

8.154 

48489.2 

1.318 

5768 . 8 

10 

2.0 

"20  000.0 

l .  53o 

11509.5 

1  .403 

"8739.5 

“ . .673 

"3516.9 

20 

2.0 

20000.0 

.  67o 

9519.3 

.574 

5964.7 

.4^3 

2348.8 

40 

2.0 

'20  0  00.0 

.533 

5  704,0 

.478 

""  2820.3 

7313  ‘ 

1607.3 

80 

2.0 

20000.0 

.323 

3447.2 

.275 

2174.1 

.18  7 

1215.7 

160 

2.0 

20000.0 

.2^0 

2565.5 

.188 

1733.3 

7119  - 

v  /  0  .  U 

5 

'4T0~ 

0  0  0  0 .0 

. 11.363 

7210.0 

1 1 .383 

721070" 

"  ”1,345 

"2364 . 1 

10 

4.0 

20009.0 

2.124 

10787.4 

1.853 

6984 . 0 

1  .278 

1923.5 

20 

4.0 

20000.0 

'  "2.504 

4316.0 

1 . 946 

2762.6 

17  065  ' 

1366.4 

.  40 

4.0 

20000.0 

1.350 

5245.4 

1.105 

2717.0 

.708 

803.0 

80 

4 . 0 

20000.0 

.898 

27 15.7 

.750 

"'  1513.0 

•  b  3  8 

494 . 1 

160 

“.0 

20000.0 

.556 

1517.9 

.390 

888.6 

.231 

460.7 

5 

8 .0 

"20000.0 

• 

*8.075 

3630.2 

*8.075 

- 3630.2 

"4.817 

— 135  7.9 

10 

8.0 

20000.0 

6.466 

3679.6 

5.133 

2740.0 

2.575 

841  .6 

20 

6 . 0 

20000 .0 

3.893 

2893.1 

3.47  1 

1 759.6 

1  .9  78 

624.2 

40 

6 . 0 

20000.0 

2.134 

1166.2 

1.818 

784.9 

1  .378 

338.6 

80 

8.0 

"200  00.0 

1.845 

1013.2 

T .  364 

602.3 

"  .97  3 

~  363.7 

160 

6.0 

20000.0 

1.186 

580.0 

1.006 

298.8 

.682 

144.8 

— 

U7 

— 

- - 

. . . . ■— 

- 

standard  error 


*#«***median  ranks****** 


*<mm*VAR  TABLES**** 

PERCENTAGE 

SUSPENSION 

sample 

6094 

40% 

0% 

SIZE 

BETA 

THETA 

BETA 

THETA 

BETA 

THETA 

BETA 

THE1  A 

5 

“5"'' 

‘  40000.0 

.680 

‘5409179.8 

. . ‘.68  0 

5409179.8  ■“ 

.  ??9 

“  52061 .7 

10 

.5 

40000.0 

.410 

1 3 1  9 1  °  7 . 6 

.35? 

424512.8 

.  )**0 

36069.2 

20 

.6 

40000.0 

.'34  b 

476201 .3 

♦  266 

130973.1 

.06  7 

24065.1 

40 

.5 

40000.0 

.161 

469644.3 

.157 

97065.7 

.092 

1117  7.1 

80 

.  5 

'40000.0 

. -.  100 

29732.7 

“.  075 

17214.4 

.087 

10053.3 

160 

•  5 

40000.0 

.067 

34843.0 

.052 

18762.3 

.0  39 

7o4b .  1 

5 

1.0 

40000.0 

3.844 

116386.8 

3.544 

1 1638R.8 

.4)7 

19065.9 

10 

1  .0 

40000.0 

.772 

6  5 2  6  3 . 3 

.600 

29591 .9 

.  28  5 

899  8.9 

20 

1.0 

40000.0 

.472 

55380.8 

.358 

J  9400 .4 

. )  8  b 

9565.2 

4  0 

1.0 

40000.0; 

.294 

10  30  17.4 

.245 

3 1 355.0 

.176 

8816. 1 

80 

1.0 

40000.0 

i  .  1  7o 

14627.2 

.146 

6813,2 

.1)2 

3830 . 2 

160 

1.0 

40000.0 

.153 

11464.8 

.125 

6546.4 

.076 

30b8. 7 

. s . 

'2.0' 

40000.0 

12. 036 

. . 28133.4 

12.033 

‘  23133.4  - 

1.591 

‘9880.9 

10 

2.0 

40000.0 

1.320 

21814.6 

1.245 

15121.0 

.718 

6315.2 

20 

2.0 

40000.0 

1.279 

16104.0 

o  6  6  b 

9642.7 

.313 

4  51-4. 9 

40 

2.0 

40000.0 

.665 

7139.4 

.620 

4695.9 

.3  17 

2417.0 

80 

2.0 

40  0  0  0.0 

.404 

‘  9906.3 

.352 

4795.1 

‘  .252 

2533.4 

160 

2.0 

40000.0 

.343 

5128.3 

.292 

3254 . 2 

.190 

1629.3 

5 

4.0 

40000.0 

“"9.517 

41023.0 

9.517 

"  41023.0 

2.203 

449  7.6 

10 

4.0 

40000.0 

1.774 

16396.6 

1.539 

1  1535.6 

1  .0**3 

3908.5 

"  20  ■" 

4.0 

40000.0 

1.641 

7861 . 1 

T  .334 

4376.2 

.831 

1516.2 

40 

4.0 

40000.0 

1.534 

4214.5 

1 .366 

2758.9 

•  686 

1210.8 

80 

4 . 0 

40000.0 

“  .766 

-  3794.4 

‘  .645 

2413.4 

. . .417 

1125.1 

160 

4.0 

40000. 0 

.523 

3333.1 

.413 

1865.0 

.276 

733.3 

''  5 

f“.'0 

40000.0 

“23.824 

.  4199.7 

23.824 

4199.7 

5.203 

2008.4 

10 

8.0 

40000.0 

8.477 

4951.3 

5.506 

3639.9 

2.053 

1692 .  ** 

20 

8.1)- 

-  40000.0 

- 32544 

-  3389.1 

"  2.470 

245  7.0 

1  .36  7 

1  259 . 5 

40 

8.0 

40000.0 

1.663 

2001 .8 

1.594 

1436.1 

1.255 

9  32.8 

80 

6.0“ 

“‘40000.0!  1.215 

1561.2 

.918 

1 054.5 

•  63  f 

570.8 

160 

*.0 

40000.0 

i  1.198 

1135.0 

.892 

667.0 

.593 

321.6 

STANDARD  ERROR 


******MEDIAN  RANKS****** 


«***VARIABLES**** 

PERCENTAGE  SUSPENSION 

SAMPL 
~  SIZE 

E 

60S 

40% 

0% 

BETA" 

THETA 

BETA 

THETA 

BETA 

THETA 

'BETA 

THETA 

"5  " 

"  .5 

“"8000  0 .0 

7.390 

29910863.2 

~  7.390 

299] 0863.2  “ 

.234 

82799.2 

10 

.6 

80000.0 

.209 

37858312.6 

.178 

4364996 , 5 

.  1  24 

4881-+. 0 

'  20 

.s- 

80000.0 

...  ...  >203 

1003922.5 

.160 

273959.9 

"  .103 

61998.9 

40  . 

.5 

80000.0 

.  142 

109270.5 

.120 

53452.4 

.069 

26439.0 

g0— 

'  80000.0 

- ,107 

. 94227.6 

.091 

. 50023.6 

.057 

243 12. 4 

1  60 

.5 

80000*0 

.  0  59 

3 73« 0.5 

.045 

22030.7 

.031 

l?27v,4 

5 

i  .o . 

i 

! 

"60000.0 

"'.866 

195172.1 

.  668 

1 95 172.1 

.8  79 

35504.5 

TO 

1.0 

60000.0 

.566 

512043.2 

•  566 

226288.4 

.  3  0  9 

28124.0 

~  '20 

1.0 

"  "80000 . 0 

.342 

'“"TO  8812.4 

“  .258 

‘  58619.0 

"‘  .209 

242 1 6.7 

\  4  0 

-w 

6  0  0  0  0 . 0 

.529 

9  4  5  ->  2 . 6 

.369 

40 1 15.2 

.  1  79 

16185.8 

1.0 

80000.0 

.249 

"42919. 1 

.194 

23175.6 

.116 

1316  7.3 

160 

1.0 

60000.0 

.202 

58800.3 

.162 

23877.5 

.  1  06 

5  8  3  U  .  1 

5 

2.0 

1 

80000.0 

28.943 

161076.5 

28.943 

161076.5 

.76  0 

1  7046.5 

10 

8.0  - 

^80  0  0  0.0 

1.276 

56147.0 

.972 

35542.1 

.499 

13678.6 

20 

2.0 

60000.0 

.  84V 

31398.3 

» b  1  7 

1 6 3 75.5 

.4]  7 

7  9..  0.4 

40 

2.0 

80000.0  ] 

.700 

30452.0 

.592 

1 6432 . 2 

.357 

6491.6 

80 

2.0 

60000.0  ! 

.510 

23486. 1 

.417 

1 2563 . 6 

•  2^D 

4  7  2  o  .  6 

160 

2.0 

60000.0 

•  366 

15841.2 

.331 

8815.5 

.  2  0  3 

2423.8 

5 

4.0 

6000 0.0 

49 .804 

40040.2 

49.804 

40040.2 

2 . 2  i  4 

“'8941  .v 

10  • 

4.0 

80000.0 

-2.202 

3fe 1 40 . 8 

1.990 

23760.5 

1  *  35U 

6d35 . 4 

20 

4.0 

80000.0 

1.43d 

1 1325.8 

1.063 

7574.4 

.598 

3411.7 

40 

4.0 

60000.0 

1 .062 

14607.6 

.652 

7936.8 

.  653 

3  0  v  3 . 0 

80 

4 . 0 

80000.0 

.  8h  1 

7495.4 

.700 

4723.8 

.4d5 

2696 . 8 

160 

!  £• 

• 

,o 

80000.0 

.626 

4288.8 

.511 

2557.0 

.323 

1708.7 

5 

«.o 

80000.0 

31.222 

13820.2 

31.222 

13820.2 

3.588 

6532.2 

10 

6.0 

80000.0 

6 . 86d 

8436.7 

6.618 

559] .1 

2.639 

3022.  7 

20 

8.0 

80000.0 

3 . 1  v  7 

9277.5 

2.240 

5683.8 

1 .3d0 

2  0  6  6 . 2 

40 

8 . 0 

80000.0 

2.102 

3953.7 

1.818 

27b 1 .8 

1.165 

1 5 1  8 . 3 

80 

8.0 

80000.0 

1.22v 

3727.6 

1.051 

2410.6 

.724 

1226.8 

160 

^.0 

80000.0 

1  .454 

2883.4 

1.142 

1326.8 

.732 

706.7 

STANDARD  ERROR 


ft  $  a  •;*  M  £  0 1  A  N  R  A  N  K  S  e  c  *  ft 

ftftftoVARIABLRS**** 

PERCENTAGE  SUSPENSION 

sample 

60*  40%  0% 

SIZE  HtTA  THE  T  A 

gpfA  THETA  BETA  THETA  BETA  THETA 

5  "75'  16000  0.0  - 

10  .  5  16  0 00  0. 0 

20  .5  160000.0 

40  .5  160000.0 

BO . .5  160000.0 

160  .6  160000.0 


5  1 . 6 

10  l.o 
20  ./  1.0 
4  0  1 . 0 

80  1.0 

160  1  . 0 


S'  2.0 
10  2.0 
20  8.0 
40  8.0 

80  8.0 
16  0  2.0 


5  4 .  o 

10  4,0 

20  '  4.0 

4  0  4 . 0 

"BO . "4.0 

160  4.0 


5  __”  6 . 0 

10  6.0 
20  "  6.0 
4  0  6.0 

80  6.0 

160  6.0 


160000.0 
160000.0 
160000.0 
i  4  o  0  0  0 . 0 
160000.0  : 
160000.0  i 


160000.0 
160000.0 
160000.0 
160000.0 
1 60  0  0 1 1 . 0 
160000.0 


160000.0 ' 
ioOOOQ.O 
160000.0  ' 
160000.0  : 
160000.0  [ 
160000.0  ; 


160000.0  i 
160000.0  ! 
1 6  0 0  0  0.0  i 
IhOOOO.O  ! 
1 600  00 . 0  I 
160000.0  i 


3.311 

« 1246581 . 1 

.47  6 

222764 ) .6 

.241 

4906541 .5 

.  1  'it 

1213650.2 

.132 

532660.6 

.0  55 

69208,6 

9  ,  h  6  9 

310432.0 

2 . 4  1  6> 

157754.9 

.597 

361878.0 

.  27)4 

953  73#  7) 

.285 

72302.2 

.l2o 

46161 . 1 

6') .  2 58  "" 

60956.4 

1.404 

95124.7 

1.211 

"  84336.1 

.611 

53747.4 

.371 

28052.7 

.  33u 

27346.7 

26.801 

38207.7 

2.526 

34804 . 1 

2 .  O  6  6 

25838 . 9 

1.045 

24249.5 

.991 

1  6  2  9  3 . 1 

.521 

100  7  H .  9 

56.569 

231 09.8 

9.407 

19404,3 

3.8  38 

18580.4 

2.677 

1 0456 , 0 

1.725 

8367.5 

1.316 

6064.5 

3.31  1 

ft  1 24858 I • 

1 

.31? 

8  h 1 395, 

6 

.181 

862046, 

.8 

.  149 

333768, 

.  1 

084 

1  3  5  9  9  d  . 

.0 

.052 

3  9  6  6  6  , 

,7 

9.469 

310+32.0 

1.341 

08974 .0 

.  349 

19} 644 .4 

.  J  82 

4  6  7  7  0,3 

.145 

41192.2 

.090 

25363.5 

68.29? 

'"  60956.4 

1.315 

59622.5 

1 .003 

44873.2 

.529 

32593.1 

'  .309 

16858,7 

.268 

15539.4 

28.801 

38207.7 

2.414 

?6?D?  .  5 

1.933 

1 7444.1 

.848 

14363 . 1 

.666 

-  7912.7 

.351 

6025.2 

56 .589 

23109.8 

3.983 

1  4  8  6  2.3 

"3.29? 

'  1  1496.2 

2,437 

5987.8 

1.379 

5343.2 

1.101 

3655.2 

.14  2 
.109 
.  1  4  3 

.090 
«  0  4  6 
.04  0 


.698 
•  3  6 6 
~  .8)2 
.  )  1*9 
.loti 
.06V 


.  6  V  0 
.  b  4  S 
.  4  4  0 
•  3  m  0 
.199 
.  1  7  2 


2.217 
.766 
.976 
.  b  v  1 
.3  9? 
.261 


6 .021 
2.43  V 
2. 18  3 
1  . ''  6  6 
•  7  9  3 
.7^2 


T  9837  6 . 8 
91535.8 
106242.0 
75348. 1 
42007.3 
27296.7 


79760.6 
S0S9S <■ v 
51266.0 
1782'.;.  1 
220)  7.9 
1 3031 . ) 


35750. 6 
23416.2 
23704.,: 
111-7.2 
7 75V. 7 
6360 . 1 


21546.) 
1  1 4  o  6 , 3 
9693.2 
7  768.2 
'4500 . 6 
32? o .  0 


7973.0 
50  3  7.5 

4380.3 

3096 . 3 
2846.0 
1  7  6  d  .  9 
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APPENDIX  IV 


RESULTS  OF  SIMULATION  STUDY  (PLOTTED) 


COMPARISON  OF  STANDARD  ERRORS 
EXPERIENCED  IN  CALCULATING  ft 
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FIGURE  17 


COMPARISON  OF  STANDARD  JURORS 

/X 

EXPERIENCED  IN  CALCULATING  & 


DISTRIBUTION  LIST 


NO.  OF 

ADDRESSEE  COPIES 


Commander 

U.S.  Army  Tank -Automotive  Command 
Warren,  Michigan  ij8090,  ATTN: 

Research,  Development  ft  Engineering  Directorate,  AMSTA-R  1 

Chief  Scientist,  AMSTA-CL  1 

Systems  Development  Division,  AMSTA-RE  1 

Transport  Vehicle  Branch,  AMSTA-REB  1 

Foreign  Intelligence  Office,  AMSTA-RI  .  1 

Technical  Data  Division,  AMSTA-RS  1 

Surface  Mobility  Division,  AMSTA-RIJ  1 

Frame,  Suspension  ft  Track  Functions,  AMSTA-R UT  1 

Propulsion  Systems  Division,  AMSTA-RG  1 

Diagnostic  Equipment  Function,  AMSTA-RGD  3 

Armor,  Materials  ft  Components  Division,  AMSTA-RK  1$ 
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